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A Theory of Ultrasonic Absorption in Unassociated Liquids 


By EE. BAUER 
Cavendish Laboratory, Cambridge 


Communicated by Sir Lawrence Bragg; MS. received 28th Fune 1948 


ABSTRACT. In many liquids the absorption of sound waves is very much larger than 
can be explained classically, in terms of viscosity and heat conduction. It is now fairly well 
established that different types of liquids behave very differently in this respect, and that 
the cause of the excess. absorption in unassociated, polyatomic liquids like benzene and 
carbon disulphide lies in the incomplete excitation, by the sound wave, of the vibrational 
degrees of freedom. It is shown that this effect leads to mean dispersion frequencies in the 
range 10° to 10!° c/s. at 300° kK. This frequency range shows that, as in the case of gases, 
the efficiency of collisions in producing excitation or de-excitation of vibrational degrees cf 
freedom is fairly low ; a theoretical treatment of this effect is outlined. Detailed consider- 
ations of temperature coefficients of absorption lead to the conclusion that all the vibrational 
levels of a liquid do not relax at the same frequency but that high excited levels relax at 
relatively low frequencies. In the case of carbon disulphide a numerical estimate of these 
individual dispersion frequencies is made. It is also shown that the variation of absorption 
with concentration in binary mixtures can be described quite simply in terms of the analysis 
developed in this paper. 


NG RO Dita LOIN 


N most liquids the absorption coefficient of sound waves, measured per unit 
length of path, varies as the square of the frequency. This frequency 
dependence agrees with that predicted assuming the absorption to be due 

to viscosity (Stokes 1845) and thermal conductivity (Kirchhoff 1868) but its 
magnitude is frequently many times larger than that calculated on these classical 
theories. It is now reasonably well established that there are at least two distinct 
mechanisms, additional to viscosity and thermal conductivity, which contribute 
to the absorption of sound waves in liquids. ‘These are the relaxation of a struc- 
tural perturbation and of some internal degrees of freedom of the molecules of 
the liquid. 

There is a definite correspondence between the chemical structure of a liquid 
and its absorption coefficient for sound waves (cf. Pinkerton 1947, 1949). The 
classical (Stokes—Kirchhoff) value of the absorption is observed in monatomic 
liquids like mercury and also in highly viscous liquids like glycerol, in which 
viscosity is the predominating factor. Water and the alcohols are typical of 
another important class of liquids which are highly associated so that the liquid 
exhibits a considerable degree of order. Here the observed absorption is up 
to three times the calculated viscous absorption and their ratio remains remarkably 
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constant over a wide range of temperature. The excess absorption in these 
liquids has been explained in terms of their relatively close-bound structure, a 
comparatively long time being required to re-establish equilibrium after the struc- 
ture has been distorted (Hall 1948). In acetic acid there is a large fall in the 
absorption and also a dispersion near 25 Mc/s. This is probably due to the 
perturbation of a chemical equilibrium, possibly the dimerization equilibrium. 
In chemically simple, unassociated polyatomic liquids like benzene, carbon 
disulphide, and chloroform the absorption coefficient has many times its classical 
value and frequently increases with the temperature. It seems likely that this 
effect is analogous to that observed in polyatomic gases, in which at room tempera- 
ture a dispersion is observed in the range 10%—10° c/s. This has been shown to 
be due to the fact that equilibrium between internal (vibrational) and external 
{translational and rotational) degrees of freedom is established relatively slowly 
(Richards 1939). Consequently, at low frequencies the sound wave is absorbed 
as a result of incomplete excitation of the vibrational degrees of freedom, while 
at higher frequencies there is a dispersion. At frequencies well above the dis- 
persion frequency the vibrational modes do not participate in the propagation 
of the sound wave. 

There have been numerous treatments of this relaxation effect for a perfect 
gas and equivalent results have been obtained by three distinct methods. 
Herzfeld and Rice (1928) introduced the relaxation process in a phenomenological 
way which does not indicate its molecular significance. ‘The analysis of Bourgin 
(1936) is based on considerations of detailed balancing and molecular transport, 
while Einstein (1920), Kneser (1931) and Rutgers (1933) developed a treatment 
which depends on thermodynamics, the equation of state of a perfect gas, and on 
detailed balance. 

Herzfeld (1941) showed that it is possible to obtain an expression for the 
absorption without making use of an equation of state. In the present paper 
a similar, but rather more detailed, thermodynamical treatment is given, with 
the aim of indicating the physical processes and assumptions involved. 


§2. THERMODYNAMICAL ANALYSIS 


In order to indicate the method of approach and to show clearly the nature 
of the basic approximations used, a general outline of the method will be given 
first. The absorption arises because the passage of the ultrasonic waves disturbs 
the statistical equilibrium of the vibrational levels, so that the effective vibrational 
specific heat depends on time. Essentially a perturbation method is used, the 
velocity of sound in the liquid being considered under conditions of disturbed 
equilibrium. A convenient way of isolating the perturbation in the expression 
for the velocity is to consider plane sound waves represented by a complex 
exponential. ‘The perturbation depends on time, so that time derivatives arise 
in the calculation: these are imaginary numbers and thus one sees that the absorp- 
tion depends essentially on the imaginary part of the expression for the velocity 
of sound. 

In transforming the expression for the velocity of sound, thermodynamic 
results are used. If there is some absorption of the wave, its propagation is not 
adiabatic or thermodynamically reversible, as is implicitly assumed when thermo- 
dynamic relations are applied. This kind of difficulty frequently arises when 
thermodynamics is applied to actual physical processes which are not quasi-static 
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and the condition for the results to be valid is that the fractional energy loss 
per cycle shall be small. Provided this condition is satisfied, the fractional 
energy loss per cycle is given by 2, where wis the amplitude absorption coefficient 
per wavelength. is given by the product of the amplitude absorption coefficient 
per unit length of path «, and the wavelength A in the liquid. For example, 
‘in benzene 24=1-6% at 10 Mc/s.; as in all liquids of this type, is proportional 
to the frequency. The same difficulty occurs in considering gases, but the 
agreement between observed and calculated absorption is quite good over the 
whole dispersive range of frequencies. 


The velocity of sound, I’, is expressed by the Euler—Laplace relation 


Tea OPISP)adiens) slam alu erates fet lar nes (1) 
where p is the pressure and p the density of the fluid. For a plane harmonic 


wave I’ is complex, and if « is the absorption coefficient per unit length of a wave 
propagated in the x-direction, then the wave motion is represented by 


exp {72av(t —x/T)} =exp {12mv(t — x/c)} exp {-—ax!} 


so that a= — mvt Ei gents (2) 


where v is the frequency, ¢ the measured velocity of sound which is equal to the 
real part of I and .4%(z) denotes the imaginary part of z. Because « is proportional 
to the square of the frequency, it is usual to consider «/v? which is independent 
of frequency. 

For an adiabatic change, the first law of thermodynamics is written in terms 
_of the internal energy E and the enthalpy H 


dE +pdsV =0; dH —Vép=0 
where E=E'+2,W,a,; H=H’' +%,W,a,, 


_a@, is the fraction of molecules in the jth excited level of energy W,, and primes 
_ refer to quantities excluding the contributions of the excited levels. All extensive 
| quantities (E, H, W;, V) are referred to one gram molecule as unit. 

It is convenient to take the change in temperature produced at any point in 
the liquid by the sound wave as independent variable, and thus equations (1) 
and (3) give 


»_ p(BH/aT)s. 
p(0E/ 0T)s 
_ Using standard thermodynamical relations this gives 


_p GOT 

aS OF TAME 

where / is the coefficient of thermal expansion, « is the isothermal compressibility, 

_ and C,, and C, are the principal specific heats. ‘The division of the thermodyna- 

mical functions into the contributions of external and internal degrees of freedom 

must be made at this stage, because it is assumed that the entropy of the whole 

assembly, S=S’ +%,a,S;, is conserved.* Thus we have 

Ee 1 (C+ 2,W;(8a;/8T)) eee (5) 
px (Ch + 2, W,(5a,/8T)) 


2 


2 


* J am indebted to Professor Herzfeld for pointing this out. 
11-2 
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The vibrational specific heat, &,W, (8a,/S7), is the only quantity contributing 
to I’ which depends on time, so that with the present formulation it is the only 
complex quantity. Taking the vibrational specific heat to be much smaller 
than C), and C’, although this is not always the case, equations (2) and (5) give 


ta ENCE o Pa Sate T Ny ‘ree (6) 


v mite 


where y is the ratio of the principal specific heats. 


83. THE REVAXING SPECIPIC HEA I 


Because of the relatively high frequency of the sound wave the specific heat, 
x,W,(8a,/5T), differs slightly from its equilibrrum value. Its value can be 
calculated most conveniently by setting up an equation of detailed balance and 
considering it under conditions of slightly perturbed equilibrium. 

First, we shall consider a single excited level, j, in equilibrium with the 
ground state which is denoted by the suffix 0. Thus a,, a) (=1—a,) are the frac- 
tions of molecules that are excited and unexcited respectively, and the equation 
of detailed balance may be written 


Od JOL= Ayko; Ghigo, eee (7) 
where the k,, are rate constants. In equilibrium, da,/dt =0, and therefore 
a; |ay =o;|Rjo = exp (— W;/RT). Rade. (8) 


The perturbation is produced by plane harmonic sound waves of frequency 
v=qw/27. Thus writing zw for 0/0t we get, from equations (7) and (8) 


§(8a,/0t) =iw3a,; =8(apkg; — ajo) =2jo3( oko; Rio — 2%) 


provided the perturbation is small. ‘Taking the temperature as independent 
variable, we obtain (Rutgers 1933) 


Lig a,| 1 + a 
wT = ERY eel 


If we disregard the excitation of harmonics of the vibration it is consistent 
to neglect exp (— W;/RT) compared with unity. ‘This gives 


da; Wi; 


xT = pall +iw/Rjol “exp (WT) eee (9) 
and if the equilibrium value of the vibrational specific heat is denoted by C;,, 
I[W,da,/5T] = —(r/y) CG [1 +0/yPF Fw. (10) 


where y = Ryg|2nr i is the dispersion frequency of the jth level. 
In ine numerical applications which are given later, the Einstein specific heat 
R(W,/RT) exp (W;/RT)[exp (W;/RT)—1]-? is used as C;. It can be seen that, 


when exp(—W;,/RT) is neglected compared with unity, the Einstein specific — 


heat.reduces to the simple form of C; given in equations (9) and (10). In the case 
of polyatomic substances there are always several vibrational levels, and equations 


corresponding to (7) and (8) can be written down for the case of / excited levels. 


‘They are 
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l 
a; =a, exp (— W;/RT); 0a,;/Ot = wedi [ages ak] -0200%. Cay 
where j=1, 2,.../ and a fraction a, of the molecules is excited math energy 


W,; k,; 1s the rate constant corresponding to an 7 transition. 
Proceeding as before, we find (/+1) linear simultaneous equations in the 
(/+1) variables 5a,/57 ; 


[( = kj) +iw]8a,/5T — % hydai8T = © ajky(W; WRT? 
tAj tAj tAj 


| 
: |. 
es hi eae (12). 
S 3a,/8T =0. 
i=0 


unknown quantities (see §5), it is not practicable to solve equations (12). — 
Because there are / excited levels, the vibrational.specific heat can be expressed 


 (Kronig 1938) as a sum of / terms of the type 


é GC, rC. 
> rece rs ——- ee eee 
j=l 1v/9; 1+2v/v, (13) 


| where the most important contribution to vy, 18 kjo/27. Results have generally been 
) expressed in terms of a mean dispersion frequency v,, defined by equation (13); 
_ from various considerations it appears very improbable that in liquids like carbon 


disulphide and benzene all vibrational levels relax at the same frequency, so that 
Vv, has no detailed physical meaning, although it does give the order of magnitude 
of the dispersive frequency range in a liquid. 
From equations (6) and (13), the absorption coefficient is given by 
(a) ae Ge 2 1 
Boze* Co vm Leu). 
Equation (9) shows that when v is very much larger than v,, the effective vibrational 
specific heat 4,W,5a,/5T becomes negligibly small, sad so primed quantities, 
which do not include the contribution of the vibrational levels, correspond to 
the case v>v 


m’* 


$4. THE MEAN DISPERSION FREQUENCY 


Equation (14) can be used to give values of the mean dispersion frequency 
for different liquids by substituting experimental values of «/v?. The results 
for a number of liquids are recorded in Table 1. 

What has actually been calculated is j,,, the absorption coefficient per wave- 
length at v=v,,, which is given by 


m 


and the mean dispersion frequency »,, is calculated from the relation 


Pane ettini LE MOY) 22 roel, ord ant hye tine (16) 
which was given by Kneser (1938) and has been used by Pinkerton (1949) to obtain 
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approximate estimates for v,,. However, both from experimental results 
(Rapuano 1947) and from a consideration of temperature coefficients of.absorption — 
(§6), as well as of transition probabilities which determine the dispersion 
frequencies y, of the. individual vibrational levels (§5) it appears that all 
levels do not in fact relax together, so that while equation (15) defines some quantity 
u* which, when substituted in (16) gives the correct v,,, yet u* isnot the absorption — 
coefiicient per wavelength at v=v,,. It follows that 2u* does not give the frac- 
tional energy loss per cycle at dispersion and thus the method used by Pinkerton. 
(1949) in estimating v,,, while interesting and good to an order of magnitude, 
cannot be regarded as supported in detail by the present considerations. 


Table1. Mean Dispersion Frequencies of some Liquids, calculated from 
Equation (16). All values refer to 300° k. 


a/v? c es C ps ried vm 
Serra (1) CA aC 6) ©) 
Carbon disulphide 10,800 (a) 1:14 (f) 1:25 4-00 0-644 Sys 
Benzene 830 (6) 1:31 (8) 153 ite? 1-300 1,200 
Carbon tetrachloride 513 (6) - 0-910 (6) 1-09 1-3 £-593 3,420 
Chloroform 380 (c) 0-987 (c) 110 7:84 0-835 2,230 
Toluene 84 (c) 1:31 (ce) 1-44. (13-0) 0-823 7,500, 
Acetone COe) - iy (ey ion 7°50 (e) 0-614 8,720 
Ethyl bromide SOO) O82 © 1-07 (8-0) 1-137 21,800 


(1) Experimental values of absorption a/v? (sec? cm! x 1071")7. 

(2) Experimental values of velocity c (cm. sec™ x 10°) (v<v). 

(3) Calculated value of c’, the velocity at frequencies v > vy, (cm. sec~ x 10°). 
(4) Vibrational specific heat C; (cal/mole deg.) (d). 

(5) Values of 2s, from equation (15). 

(6) Mean dispersion frequency vy, (Mc/s.). 


(a) Kneser (1938). (d) Herzberg (1945). 
(6) Pellam and Galt (1946). (e) Schumann and Aston (1938). 
(c) Willard (1941). (f) Biquard (1935). 


+ The value of a/v? due to viscosity and heat conduction has been subtracted. 


It may be convenient to give here the formula for the dependence of the 
velocity on frequency, which is of course the real “ dispersion” formula. From 
equations (4) and (8) we obtain 


uC; y'-—l 
| is ial j 
C=C E itG)ae © ] sei ee (17) 


since 1,W,da,;/5T<C, C,. This equation shows that lim c,=c’, as has already 


been inferred from physical considerations. 


§5. THEORETICAL ESTIMATES OF THE DISPERSION FREQUENCY 


From equations (9) and (10) it can be seen that to calculate y, the essential 
step is the calculation of k,o, the rate at which a molecule in the jth excited vibra- 
tional level falls back into the ground state. There is no reason to believe that 
there exists a potential barrier large compared with RT to prevent de-excitation 
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| (Eucken and Becker 1934), so that the rate constant kj) is given by the product 


| of the collision frequency y,,, and the transition probability p,). For a gas, the 


» collision frequency can be calculated from the mean free path and mean velocity 


_ ofthe molecules; at room temperature and atmospheric pressure v,4~ 5-109 sec}. 
For liquids, vj ~RT/h (Glasstone, Laidler and Eyring 1941).which is of the 


+ order 102 sec-! at 300°K. The observed dispersion frequencies for gases lie 


| in the range 10°-10%c/s., and predicted values for liquids are of the order 
10°-10'° c/s. Thus the transition probabilities lie in the ranges 10-® to 10-3 
for gases and 10-* to 10-? for liquids. 

The low value of this transition probability in the case of a gas has been 
explained by Zener (1931), who showed that it is a consequence of the quantal 
treatment of collisions between an atom and a vibrating molecule that inelastic 
collisions are infrequent. Later, Landau and Teller (1936) derived an approxi- 
mate classical formula for the dispersion frequency, but their result is not applicable 
at room temperature where the spacing of the vibrational levels is relatively 
large compared with RT. 

In the case of liquids it is easy to modify Zener’s treatment by adding a constant 
attractive potential Up, representing the cohesion, to the Lennard-Jones repulsion 
between molecules, Ey exp [—«(r—7,)] which of course depends very strongly 
on r, the distance apart of the molecules. ‘Then, if ¢,; is the energy of the jth level 
and K is the kinetic energy of the colliding particles before the impact, the transi- 
} tion probability is given by 

ar Ke (OF Ke Us, eae 
we 16 ponies “(eG Upts; 
where all quantities are measured in atomic units, and w and p,, are reduced 
masses of the system as a whole and the vibrating molecule respectively. An 
outline of the derivation of (18) is given in the appendix where order of magnitude 
values for the various quantities are indicated. .They give values of pj) which 
_ are of the correct order (10-° to 10-?) and in view of the idealizations and approxima- 
tions of the analysis no detailed deductions can usefully be made; in particular, 
the difference in the behaviour of different substances cannot be explained 
because no account is taken of the structure of the colliding particles or of dipole 
moments which produce an increase in the transition probabilities (Metter 1937). 
From equation (18) it can be seen that p,. decreases with increasing energy of 
excitation e, and, as both K and Upare of the order 3R7/2, the transition probability 
increases with j increasing temperature. 

Eucken and Becker (1934) found that the transition probability in gaseous 
mixtures increases as the third or fourth power of the temperature. In an 
attempt to explain the inefficiency of collisions they showed that if the efficiency 
of a collision is inversely proportional to the 2th power of the closest distance 
of approach of the molecules, then this transition probability is proportional 
to JT”. Exactly the same temperature variation results if one supposes that the 
efficiency of a collision is proportional to the mth power of the mutual kinetic 
energy of the particles before impact. In the case of liquids this picture again 
has a qualitative significance, but because of the interaction represented by Us, 
the value of n will be rather smaller, of the order 1 to 2. 
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§6. THE TEMPERATURE COEFFICIENT OF ABSORPTION 


We write equation (14) in the form 


ie a v i 
const nea ie tea ae (14a) 
Vv ; c Vin 


where the constant, zy/y'C), is approximately independent of temperature, 
while by a standard thermodynamic result (y’—1) is proportional to 7. Now 
it is convenient to write 0 log C,/2T =D,/T where D, ~ W,/RT—2 for W,>RT 
and 0 log £C,/dT = D/T. It has beea seen that the transition probability increases 
with temperature, approximately as 7”, where 7 is of order 1 to 2. The collision 
frequency, RT/h, is proportional to 7, and logarithmic differentiation of equation 
(14 a) gives 

05 re Cen eee 

wT re EA 


By consideration of Table 2 we see that, for this equation to be satisfied, » would 


[Deis +t si aeneeee (19) 


Table 2. ‘Temperature Coefficients of Absorption 


Ome derma 10> deo= (D—D) err D=T9 log LCS Or Ni 


Substance (1) (2) (3) (4) i 

3 a 

Carbon disulphide (10) —(3-5) +1:95 0-874 
Benzene 11 (a) — 3:6 (a) +2-22 1-820 

Carbon tetrachloride (0) —(3°5) —1-05 0-623 a 

Chloroform (10) — 3-4 (b) +1-98 0-666 ; 

Acetone (10) — 3-8 (b) +1-86 1-840 : 

Ethyl bromide 10 (a) — 3:6 (a) +1-:92 (1:0) ; 


(1) Measured values of temperature coefficients of absorption (10% deg-?). 
(2) Measured values of temperature coefficients of velocity (10% deg—?). 
(3) (D—n)ere calculated from equation (19). 


(4) D=T@log UC,/OT (see § 6). 
(a) Pellam and Galt (1946). (6) Willard (1941). 


* Experimental values for temperature coefficients of absorption are rather incomplete, 
but for unassociated liquids the values are generally near (+)10~? deg-! ; carbon tetra- 
chloride is an exception, as here the absorption seems to vary very little with temperature 
(Bazulin 1937). 


have to be much less than one or negative in many cases. These values cannot 
be reconciled with equation (18), or with simple physical considerations like those 
of Eucken and Becker. 

The difficulty can be overcome if we suppose that the individual levels relax 
one by one, a high excited level having a low probability of de-excitation and thus 
a low dispersion frequency. Because a high excited level has a low specific 
heat but a large temperature coefficient of specific heat, the change in weight 
of the several levels arising from the replacement of 


log (UC;/,)//OT by dflog (= C;/y)]/aT 


would result in a larger value of (D —7),,.. 


Unfortunately there are no data available on the values of v,, but there are 


several independent results which show that in certain cases different vibrational 
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© levels relax at different frequencies. Acetaldehyde vapour (Alexander and 
_ Lambert 1942) shows three distinct dispersion regions, which the authors ascribe 
_ to the relative complexity of the molecule. 


The behaviour of liquid carbon disulphide is of special interest, because it 


is the most highly absorbing simple liquid which has been investigated. A 


study of the absorption coefficients measured at different frequencies shows that 
multiple dispersions must be assumed to explain the results which are given in 


Table 3A . 
In Table 3B these results are interpreted in the light of the previous discussion 


and the dispersion frequencies of the three vibrational levels are calculated. 


Table 3. Multiple Dispersion in Carbon Disulphide 


A: Experimental values of the absorption 
my" (cm =i sec? x 10-1”) 10,750+ 1,000 7,400 1,400 
Frequency (Mc/s.) il 6:5 75,105 


Observer Claeys, Errera, Sack Willard (1941) Rapuano (1947) 
quoted from Kneser (1938) 


B: Interpretation of these results 


Vibrational Vibrational Velocity 
level specific Cc es 2h; A(a/v?) Vj 
frequency heat, C; (v<v;) (v>v,)  [cf.eq.(15)] (cm=tsec? (Mc/s) 
(cm!) _(cal/mole deg.) (cm. sec! (cm. sec~! Sel Om) 
K 1102) x 10°) 
1523 0-0724 1:14 1-14 0:00705 3,350 1:85 
656,5 0-930 1:14 iletS7 0-106 6,000 SiS 
397 3-00 tail yy) 1-245 0-494 1,400 305 
(doublet) 


It is now possible to evaluate 
evap 
aa EEG 1 
XC; /vy3 OT; 


for carbon disulphide on the basis of Table 3B. The value obtained is 2:76, 
so that if is taken equal to one, (1/«)(0«/0T)=9-10- deg". 


§7. BINARY MIXTURES 


It is well known that in the case of gases small traces of impurities may have 
an extremely large effect on the absorption (Eucken and Becker 1934). A similar 
effect has been found in a number of liquids like carbon disulphide and benzene 
where the absorption falls very rapidly when a small quantity of a less absorbing 
liquid is added to one of high absorption. Pinkerton (1949) has shown that this 
effect can be explained on a very simple relaxation picture. It is shown here 
how the same basic ideas, when interpreted in terms of transition probabilities, 
yield results which are in good agreement with experiment. 

Let A be the strongly absorbing liquid and B the other one. Then A has a 
much smaller dispersion frequency than B, so that an A molecule, once excited, 
has a much smaller chance of de-excitation than an excited B molecule. For 
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simplicity in the analysis it will be assumed that only binary collisions need be 
considered and that collisions between excited molecules may be neglected. 
Then if A and A* represent an unexcited and an excited A molecule respectively, 
and similarly for B, B*, then an (A*A) collision is much less efficient than a (B*B). 
collision because vis much smaller than v8. In a mixture of molecules of both 
species, A and B, there will, in addition to (A*A) and (B*B) collisions also be 
(A*B) and (B*A) impacts which have a rather high de-excitation efficiency 
because it is in general impossible for a quantum of vibrational energy to be 
transferred from one type of molecule to the other without some loss of energy. 
Now, following Pinkerton, we make the crucial assumption that (A*B) and 
(B*A) collisions are as efficient in producing excitation and de-excitation as. 
(B*B) collisions. Thus, as the concentration of B molecules increases, the net 
efficiency of all collisions tends rapidly to the value corresponding to the liquid B,. 
so that the absorption falls sharply. 

In the analysis it is convenient to suppose that each molecular species has. 
only a single excited level of energy W, and W, for A and B respectively. The 
generalization for a number of levels proceeds as in §3 for the case of a single: 
liquid, and the same approximation is made. 

Consider 1 gm. mol. of the mixture; a fraction a is made up of A molecules,. 
there being mole fractions a,, a9(=a—a,) of excited and unexcited molecules, 
and similarly mole fractions b,, 6) for the liquid B, where a+6,+6,=1. Intro-. 
ducing k, g, /andh respectively as (AA),(AB),(BA) and (BB) transition probabilities, 
we can write two equations of detailed balance: 


Ca 
a = (Ro Go + £0190) a0 — (R10% +1090) = $014 — $10% 


* Ss 
7 = (L014 + o140)b0 — (L104 + 21090) 01 = foro — $1051 


Proceeding as in §3, one arrives at the result 
0a,/OT =(Wy/RT*)a[1 +ie/dyo-texp(-W,/RT) | (21) 
0b, OT =(Wy/RT*)(1 —a)[1 + ieo/thyo}-t exp (— Wy/RT] in a ; 


and if C,, Cp, denote. the vibrational specific heats of A and B, equation (13) 


becomes 
8) -m[BG DALES SA) ne 


Now two approximations are introduced: 


(i) The crucial assumption that, to a first approximation, (B*B), (B*A),. | 


(A*B) collisions are all equally efficient in producing an inelastic transition, 


z times (z~10) more so than (A*A) collisions. In mathematical form this is. 


written 3k;;=g;;, =], =h;; Me 70, <1); 


a 


1 
(ii) The quantity 27? aly’ -1)a changes linearly with concentration a. 
p 


This is not a critical Ricoh because the total variation in this quantity is. 


not large. It will be written M,(1 —fa). 
Then equation (22) becomes 


(a/»*)a=(Mp/g10)(1 —fa)[(1 — a) Cy, + a2C/{2 — (2 —1)}] 
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which can conveniently be compared with experiment in the form 

CaN az C 

— =(1- 1—a)+ ——_ |. a, 

sf =(1-fa)| (1a) + — | (23) 
Results for a benzene-toluene mixture are given in the Figure. If A specifies. 


10 (6) 


O/O Toluene 
nD 
ee N 


SS 


2 4 
Ca 
a 
0 
0 02 04 06 08 10 
Mole -fraction of Benzene 


Ratio of the absorption coefficient of a mixture of benzene and toluene to the absorption 
coefficient of toluene as a function of the mole fraction of benzene. 


Full line : theoretical curve (equation (23)). 
Circles : experimental points (Grobe 1938). 


benzene and B toluene then, using the values given in Table 1, f and z are given: 

by f=1— 2am —0-743 and zg=v8/vi=6-28 since C,/C,z=0-904 and 
“BR Cr os 

a, /% ,=9-80. Pinkerton’s (1949) result is recovered by putting f=0 and 


C, =Cg so that z=a,/ap. 


§8. DISCUSSION 


Since the present approach to the problem of vibrational relaxation is based 
on rather general considerations regarding the perturbation of an equilibrium, 
the same general method can be applied to a number of different aspects of the 
problem of incomplete excitation of some internal degree of freedom of the mole-. 
cules. In particular, the application to the perturbation of a chemical equilibrium 
by the sound wave is of interest. Such an effect may be the cause of the dispersion 
observed in acetic acid by Spakovskij (1938) and others. Pinkerton (1948), 
using his own and Bazulin’s (1936) results for the absorption of sound in acetic 
acid near 2 Mc/s., has estimated by a graphical method that the fall in «/v? in 
this frequency range is 12,000 x 10-1? cm™! sec? and that the dispersion frequency 
is 2 Mc/s. at 17°c. From these data a value of the relaxing specific heat can be: 
calculated from equation (14): its value is 1-4 cal/mole degree. ‘This is very 
large and might correspond to a hindered rotation or similar process. 


§9. CONCLUSIONS 


In the case of simple, unassociated polyatomic liquids like benzene, incomplete: 
excitation of the vibrational degrees of freedom by the sound wave is shown to. 
lead to a dispersion, with mean dispersion frequencies in the range 10° to 101° c/s. 
at room temperature. It does not seem probable that all vibrational levels of 
a liquid relax together : it is possible that some high excited levels have dispersion 
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frequencies below 107 c/s. which have not so far been detected. Carbon disulphide 
may well be a case in point, for consideration of the available experimental data 
indicates that its three vibrational levels may have dispersion frequencies of 
1-85, 15-5 and 305 Mc/s.. 
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A‘P PEN DIX 
Calculation of Transition Probabilities 


A brief outline of Zener’s (1931) treatment is presented here, with a modifica- 
tion to account for the cohesive forces in a liquid. ‘To facilitate reference to the 
original paper a similar notation is employed and atomic units are used. 

Zener considers the one dimensional collision of an atom A of mass WZ, with 
a vibrating molecule B-C, the two atoms of this having masses M, and Mg. 
By a suitable choice of coordinates the motion of the centre of mass becomes 
separable and is eliminated, and the wave equation of the system is written 

2 2 

| 7 F — V(x)— U(x, +E eat eke MES Peel (24) 

Here’ is the coordinate describing the position of A relative to the centre of mass 
BC, x is the distance BC between the vibrating atoms, pu,, =MyM,/(My+™M,) 
and w= M,(M, + M,.)/(M, +My+M_,) are the reduced masses of the vibrator 
BC and the whole system respectively; E is the total energy of the system, V(x) 
the potential energy of the vibrator, and U(x, 7) the mutual potential energy of 
atom A and molecule BC. The modification which has been made to Zener’s 
treatment is that a constant attractive term — U) is added to the repulsive term 
V(x, r) which depends strongly on distance. Uy must be present, and rather 
larger than 3R7/2, for the liquid to maintain its volume. 


Ulein=Ue a ere a ns eee bee (25) 


The vibrational motion of (BC) in its 7th state is described by the wave 


equation 
[1 /Mm)(@?/dx*) — V(x) + Bi]e(x)=0 sa (26) 


‘where $,x) is the normalized wave function corresponding to the energy £,. 
Now, define the momentum, p,, by 


Pi = [mE ss Uy is E;)] t DACA SA (27) 
We require a solution of (27) having the asymptotic form 
Pit + eo Dir Ee Per 


pi ae rar a p(x). eon) sions (28) 


"This solution of course consists of an incident wave which is scattered, in part 
inelastically: the probability of a change 7—>f in vibrational quantum number 
as | Vig?» 


V(1> 0) = 


, where V;,(r) = [d,(x) V(x, r)4j(ex) de. 


A Theory of Ultrasonic Absorption in Unassociated Liquids 153, 


To obtain Yig> We must find a solution for (24) of the form (28). The 
wave function ‘is expanded as Y= LOUr) Pel), the O,(r) being undetermined. 


functions. ‘This is substituted in eehation (27) and after multiplying by ¢,(x); 
and integrating over x, the orthonormality of the ¢; leads to 


fe — HV + | Ol) = O)V fe) vee (29) 


j#k 


Born’s collision method is now used to obtain approximate solutions of (29) ;. 


| first the right hand side of (29) is put equal to zero and the zero-order solution, 


((r), corresponding to elastic scattering is found. Then Q(r) is substituted 


on the right hand side and the equation solved again, to give the first order 


approximation, O{”(7) which is the one used. 

The remainder of Zener’s treatment is purely mathematical. In order to 
obtain a relatively simple set of solutions of equation (29), he replaces the potential 
Vy.=E,exp {—«(r—79)} by an inverse square law which fits the exponential 
function in the critical region where r~r). This enables the solution of (29) 


| to be found in terms of Bessel functions and finally gives the transition probability 


for de-excitation of the first excited state of the jth level, of energy e,, 


2 
ga 7 oH Eo (Pi\VuBo- Me 
o=|AbP= 767 2a (2) Ae niente (30) 


If one writes EH —E,= K, the relative kinetic energy of the atom A relative to the 
molecule BC before impact, equation (27) gives 


wT pp Eo feK+ Oy ae 
[Kae Uy+<; : 


Substitution of numerical values shows that this expression gives results 
which are of the correct order of magnitude. ‘Taking A, B and C to be each of 
molecular weight 50, Ey=RT, K=U)=3RT/2, ¢, =3RT at 300°K. (~625 cm) 
then for « equal to 2 and 3 in atomic units the corresponding values * of pj are 
6:04 x 10-3 and 3-12 x 10-. 


Note added in proof (8th JFanuary 1949). The treatment of transition 
probabilities which is outlined in the appendix cannot be regarded as completely 
satisfactory in its application to liquids, as it is but a modification of the 
corresponding analysis for a gas. It would be better to use a quasi-crystalline 
model of a liquid, and work on this is proceeding. 


Pio i 16 Pm €j 
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Second Sound and the Behaviour of Helium II 
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ABSTRACT. In this paper some theoretical consequences of the existence of second 
sound in helium II are discussed. A method is suggested for setting up a quantitative 
phenomenological description of the thermodynamical behaviour of helium II. It is known 
that expressions for the velocities of first and second sound may be obtained by using quite 
general models. ‘The thermodynamical behaviour may be described in terms of the 
velocities of these two types of sound, so that it becomes possible, in principle, to obtain 
explicit values of the unknowns by combining the two sets of equations. On the basis of 
these ideas a discussion is givén of the coefficient of thermal expansion of helium, the 
ordinary thermal conductivity, and the attenuation of second sound. 

The paper contains a discussion of the modifications necessary to the Debye theory of 
-specific heats when the sound velocities are temperature-dependent. 


§1. PREVIOUS THEORIES OF THE THERMODYNAMICAL 
BEHAVIOUR OF HELIUM II 


anomalous properties. At 2:19°k., the A-point, a phase change of the second 
order takes place, and below this temperature the liquid, known as helium II, 
possesses properties rather similar in nature to those of superconductors. Details 
of these properties may be found in Keesom’s book on helium (Keesom 1942). 
A particular property with which we shall be concerned in this paper is that below 


[ is well known that liquid helium at very low temperatures possesses 


* Now at the Royal Society Mond Laboratory, Cambridge 
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>the d-point the values of the entropy and the specific heat of the liquid are very 


»much larger than the values to be expected from Debye’s theory (Debye 1912), 


>) calculated from the known values for the velocity of ordinary sound in the liquid. 


di Theories have been suggested by several authors to account for these properties ; 


A may mention in particular those of London (1938, 1939), Tisza (1938, 1940, 
11947), Landau (1941), and Green (1948). 


___ The ideas of London and Tisza are based on the concept of a Bose-Einstein 


«condensation, a phenomenon predicted theoretically for an ideal gas (London 1938, 
641939, Fowler and Jones 1938, Dingle 1949). According to this theory, below a 
1 certain very low temperature a fraction of the available atoms become condensed 
J) into the lowest available state. These authors believe that many of the peculiar 


) liquid helium. However, it must be stressed that the theory is only a qualitative 


ie of this ideal Bose-Einstein gas may be taken over into the theory of 


) 
: 


one, since the values of the entropy as a function of temperature given by ideal gas 


‘| theory are too inaccurate to be of value, the entropy being given as proportional to 


” T?° as opposed to the experimental relation Soc T*5. In order to obtain reason- 
qd able agreement with experiment, the model must be altered in an empirical 


| manner. 
) Landau’s general theory of helium (Landau 1941) is based on his quantum 
1 hydrodynamics, which is an attempt in quantum theory to avoid the difficulties of a 


4 complete microscopic theory of fluids in a similar way to that in which they are 


/ avoided by classical hydrodynamics. Landau concludes that only discrete 


» values of the curl of the velocity are permitted, and that there will be a gap in the 


| energy spectrum associated with these ‘“roton” excitations. The entropy will 


) be contributed by ordinary Debye waves (phonons) and by these rotons. The 


| phonon contributions are calculated in the usual manner. In order to derive the 
roton contributions, Landau assumes that only the first gap in the roton energy 
/ spectrum is of importance, and the excited levels are replaced by a continuum for 


i which the energy « is given by «-=A+p)?/2u, where A is the energy gap, p the 
? momentum, and p the effective mass of a roton. ‘Thus in Landau’s theory we 
{ have still two arbitrary parameters A and », which must be fixed by appeal to 
experiment.* - 


The very recent work of Green (1948) is based on the general kinetic theory of 
| liquids developed by himself and Born (Born and Green 1946, 1947a, b, 
Green 1947, 1948). On quantization it appears that pressure defined thermo- 
| dynamically is no longer identical with pressure defined by kinetic theory. Green 
believes that many of the peculiar properties of liquid helium may be explained 


as due to this difference. Unfortunately, the detailed theory proves too compli- 


cated for direct application, and Green’s results (e.g. the magnitude of this 
pressure difference) have in fact tc be derived from experiment rather than from 
theory. 


§2, QUALITATIVE DISCUSSION OF THE THEORY PROPOSED 
The London-Tisza model, based on the concept of Bose-Einstein condensation, 
and the Landau theory, based on quantum hydrodynamics, agree both between 
themselves and with experiment in so far as the helium IT is considered rather as a 
mixture of normal fluid and superfluid. On the London-Tisza model, the super- 
fluid is supposed to consist of atoms which have piled up in the ground state, due to 


* In a recent variant of the theory there are three parameters (Landau 1947). 
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some kind of Bose-Einstein condensation, and the normal fluid to consist of the — 
remaining atoms, i.e. those in the excited states. On the Landau theory, the — 


superfluid is likewise made up of atoms in the ground state, whilst the normal fluid 


is described by excitations of Debye waves (phonons) and excitations of vortices — 


(rotons). In both theories the entropy of the superfluid is of order zero, since 
only atoms in the ground state are involved. Tisza and Landau have shown 
independently that in a fluid of this nature we should expect two types of sound: 
ordinary “first”? sound caused by the movement of the normal and superfluid 
together, and a new type, “second” sound, caused by relative movement of 
normal and superfluid. In second sound, owing to the relative movement of 
normal fluid, which carries entropy, and superfluid, which does not, we obtain a 
state of the liquid such that alternate layers contain an excess or a deficit of entropy. 
It follows that these layers are respectively at an effectively higher and iower 
temperature. ‘Ihe most important difference between the two types of sound is 
that the dominating feature in first sound is fluctuation in pressure, whilst in 
second sound it is fluctuation in temperature. Similar conclusions have also been 
reached by Green (1948). 

On two-fluid types of model, the velocity u, of second sound is given by 
us = TS*p,/Cp,, where T is the temperature, S the entropy per gramme, C the 
specific heat per gramme and p,/p, the ratio of the densities of superfluid and 
normal fluid. In previous theories, more specialized models, or empirical 
formulae, have been used to calculate the values of S, C and p,/p,, and hence of 
u,asafunction of temperature. Inthis paper we suggest a method by which these 
special models may be avoided. ‘The argument is that a quantitative thermo- 
dynamical description of helium II may be given in terms of the velocities of the 
two kinds of sound by employing an extension of Debye’s theory of specific heat: 
this gives us the values of S, Cand p,/p, in terms of u,, assuming that these two types 
of sound provide the dominant contributions to the thermal energy. In principle 
we have only to combine these equations with that giving w, in terms of the thermo- 
dynamical quantities to enable us to determine all the unknowns. In practice it 
will be much easier to show that the experimental values of uw, are consistent with 
both sets of equations. 

We shall encounter three main difficulties. Firstly, we know that there are 
3.N degrees of freedom per unit volume, where N is the number of atoms per unit 
volume. ‘Thus there are 3N normal independent modes of vibration per unit 
volume : our problem is that of allocating these between first and second sound. 
A similar difficulty arises in the usual Debye theory, namely the question of the 
allocation of the available modes between longitudinal and transverse waves. 
‘Here some authors have made the arbitrary assumption that the frequency 
limits of the two types of waves may be taken to be the same; Brillouin (1936) and 
Lucas (1937) have shown that the assumption of equal wavelength limits is 
probably preferable. In our case this difficulty is not serious, as we shall find that 
second sound contributions are of overwhelming importance in the temperature 
region so far investigated. 

Secondly, the usual Debye theory giving the thermodynamical behaviour in 
terms of the velocities of sound is only valid if these velocities are independent of 
temperature, the substance being at constant volume. This condition is roughly 
satisfied for first sound, but not for second sound, particularly near the A-point. 
We shall discuss this problem in §§8 and 9 of this paper. 
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_ The third difficulty is that of determining how to take into account the absorp- 
‘ie and dispersion of sound waves, as it is known that these effects become 
Beportant as the A-point is approached. 


Soy EE DEBYE THEORY 


Debye (1912) considers that in solids (and liquids) waves of sound will be: 
| present continually, the frequencies represented being those up to a certain limit 
j dictated by the atomic structure of the material. This upper limit may be 


‘calculated from the fact that a real solid (or liquid) has only a finite number of 


+ degrees of freedom, so that we must cut off the upper end of the frequency spec- 
trum to make the theory correspond. Physically, this corresponds to the idea 
| that if the wavelength of the supposed Debye waves were shorter than the inter- 
} atomic distance, the vibrations of shorter wavelength would not be independent of 
| those of longer wavelength, since it is always possible to represent displacements of 
. a finite number of particles equally well by an infinite number of vibrations of 
) wavelength shorter than the interatomic distance. 

. For a single sound velocity c, independent of temperature, the Debye theory 


| yields the following values for the thermodynamical quantities in which we are 
interested (see §8): 


y 
: E T\3 (if x3 dx 
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| where here C is the specific heat (at constant volume) per molecular -volume. 
'@ is the Debye characteristic temperature given by the relation 


RUAN Actes (4) 
At sufficiently low temperatures, T7<6, we may use the simpler expressions 
GT = IR(T IO) aAf15, VANEOMEY eS PRES C8 (5) 
SES VAR IUD ici te et ee AN (6) 
De OIRO Sato id ot llc OR Bile ose (7) 


St: THE AP Page Ain LON Ok 2 hE sD rE BYP. “THEORY 


The formal theory quoted in the last section refers to the case of a single type 
» of sound wave of constant velocity c. This differs from our case in two respects: 

(a) It is an experimental fact that with decrease in temperature below the 
, A-point the velocity of first sound rises steadily with temperature,* and that of 
second sound rises rapidly from zero at the A-point to a value which may be con- 
| sidered as roughly constant ¢ between 2-0°x. and 1-03"k., the lowest temperature 
for which it has been measured, For a preliminary discussion we shall apply the 
results of the Debye theory as given above, replacing the temperature-dependent 


* Since the coefficient of expansion of helium II is both very small and negative, this holds for 


- constant volume as well as for the measured results at constant pressure. 
t More accurately, there is a very flat maximum in the curve of velocity against temperature at 


about 1-6° K. and a minimum at 1:12°K 
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sound velocities by mean values. In §§8 and 9 of this paper a detailed treatment 
of this problem will be given; this yields information on the possible shape of the 
curve of second sound velocity against temperature. 

(b) The results of the last section have been given for the case of a single type 
of sound wave of velocity c. In our case we have two types of sound: first sound 
of velocity u,, and second sound of velocity u,. The question to be decided is that 
of allocating the 3N modes of vibration available. We may discuss the following 
simplifying assumptions : (i) The frequency limits are taken to be the same, in 
which case we may introduce a mean velocity c defined by the relation 
1/c3 =1/u} + 1/u3, and then apply the theory given above without further alteration. 
We are assuming here that only longitudinal waves may pass through the liquid. 
(ii) The wavelength limits are taken to be the same, in which case we must take the 
frequency limits as proportional to the respective velocities. To preserve the 
total number of normal modes at 3N, we must take the frequency limits as 
vy =(9N/87)*u, and v,=(9N/87)*u,. We then consider the contributions from 
the two types of sound separately. In the known temperature region, where 
Uy ~U,/13, the effects of first sound may be neglected, so that assumptions (i) and 
(11) give the same result. 

There is also a third possibility, which arises in the following way. In first 
sound our frequency limits are really connected with the point at which the concept 
of hydrodynamic flow becomes meaningless; this limit will be determined by the 
condition that the wavelength becomes of the order of magnitude of the interatomic 
distance, In second sound our frequency limit is determined by the breakdown 
of the concept of relative hydrodynamic flow between the superfluid and the 
normal fluid. This breakdown will occur when the wavelength is of the order of 
magnitude of the mean distance between superfluid atoms if these are in the 
minority, or between normal atoms if these are in the minority. We may then 
show that the characteristic temperature must be proportional to (p,/p)* or to 
(p,/p)', whichever is the smaller at the temperature considered. An objection 
to this argument must be raised. When we speak of anormal or a superfluid atom, 
we must remember that the chosen atom will remain “normal”? or “superfluid” 
for a time of the order of the relaxation time between different energy states. If 
the relaxation time involved were much shorter than the reciprocal of the frequency 
of sound concerned, the avove argument breaks down, as superfluid and normal 
atoms constantly interchange, so that we again become interested in the much more 
slowly varying mean distance between atoms, regardless of the energy levels which 
they occupy. 


§5. COMPARISON WITH EXPERIMENT 


A formally complete description of the thermodynamical behaviour of 
helium II should be given by the combination of the following equations. (The 
effects of first sound will be neglected, a procedure justified in the temperature 
region so far investigated.) 


(1) The equation u} = T'S%p,/Cp, =(1—p/p,)/d(1/S)/dT. 
(2) Equation (2) of §3, or more accurately equation (11) of §9, giving S in 
terms of wp. 


(3) The equation p,/p=4E/3u5 (or more accurately equation (14) of §9), 
an extension of an equation due to Landau (1941). 
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We have now sufficient equations to eliminate the unknowns—at least in 
principle. However, as has been remarked previously, it is much simpler to take 
the experimental values of u, or, alternatively, trial values, and check that they are 
consistent with these three equations and with the experimental values of the 
entropy. 

First, it is instructive to calculate the entropy taking only first sound effects 
into account—1.e. to find the Debye entropy calculated according to previous 
theories. It is known that the velocity of first sound in helium II is about 240 
metres per second, and that the molecular volume is 27-6cm?. From these we 
deduce * the value 6,~ 29°. from the relation 0, = hy,,/Rk = hu,(9N/47)'. This 
value of the Debye characteristic temperature is so much higher than the temper- 
atures in which we are interested that we may use the simple formulae valid in the 
T? region. ‘The value of the entropy derived is 1-8 x 10-* T?cal/gm., which is to be 
compared with the experimental relation S=0-405 x (T/T )>*cal/gm., Ty being 
the temperature of the A-point. ‘Thus near the A-point the experimental value of 
the specific heat (7'd.S/dT) is roughly 40 times the value calculated from the usual 
Debye theory, taking only first sound into account. 

Let us now introduce second sound. At first sight it might appear that since 
C~1/c? (in the T* region) and u, ~ u,/13, second sound would contribute a specific 
heat about two thousand times as large as that of first sound. This is not so, 
because the introduction of second sound effects causes us to pass right out of the 
T? region for all except the very lowest temperatures. ‘This may be seen from the 
relation 0c, showing that 0. ~6,/13~2-2°K., so that 6, is of the same order of 
magnitude as the temperatures we wish to consider. For instance, to take a very 


-O/T 
simple case, the energy is proportional to the integral x? dx/(e~—1) When 
0 


6/T>1 (i.e. in the T? region), this has the value 6-45, but when 6/T=1, say, it 
takes the value 0:224. Thus instead of second sound effects giving an energy 
approximately 2,000 times that previously calculated (assuming only first sound), 
it provides a contribution approximately 70 times as large. ‘The other thermo- 
dynamical quantities are cut down in a similar way. 

An experimental fact that is particularly interesting is that the specific heat 
just above the A-point is unduly high, so much so that the slope of the curve of 
specific heat against temperature is negative. This may well imply that it is 
possible to transmit hypersonic waves of second sound at temperatures just above 
the A-point ; it may also be possible to transmit waves of attainable frequencies. 


VOL Divo RIVA EXPAN STON OP THE LIU 1 


Generally speaking, solids and liquids have a positive coefficient of thermal 
expansion. We wish to explain why the coefficient of thermal expansion of 
helium II is very small and negative in sign, whilst that of helium I is of normal 
magnitude and of normal (positive) sign. First we must find a quantitative 
relation between the coefficient of thermal expansion, «, say, and 0, the character- 
istic temperature. For the present we shall make the usual assumption that @ 

* Tisza (1947) has given the value 18° K., but we are unable to see how he arrives at this result. 
We also cannot agree with his criticism of Pickard and Simon’s (1939) experimental value of 0, since 
the formulae for the specific heat, etc., given in terms of # are completely independent of theoretical 


speculations concerning the existence or non-existence of transverse waves in liquids, these all being 
collected in the connection between 6 and the velocity of sound. 


I2-2 
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depends only on volume—i.e. 9=0(V). Referring back to (2), we see that the 
entropy is a function of T/@(V). Thus S=S(T/0(V)), so that 


po aoe eo 7) --">(F) 
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where S’ is the differential coefficient of S with respect to the variable 7/6. 
Eliminating S’, we have 


as) __ (1) (36) (as 
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Applying the Maxwell thermodynamic relationship (0S/0V),=(0p/0T), and the 
relation C= 7(0S/0T),, for the specific heat at constant volume, we obtain 


Ss he ating Sy (6oas8l Wire ad 
(5), ~ GNOV fae Ve 


where y=(0In6/dIn V), and V is the molecular volume. Now since 


ni (OP op 
dp = (sh),@T+ (sf),,a”. 
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where « is the coefficient of thermal expansion and y is the isothermal compressi- 
bility. Thus «= —yCy/V. 

To sufficient accuracy we may treat y and V as independent of temperature. 
It is worth noting that since « is proportional to C, the coefficient of expansion will 
have a peak value at the A-point, and will fall to zero at very low temperatures. 

The sign of the coefficient of expansion will be the same as the sign of —y, i.e. 
the sign of (—dc/dV),, where c is the velocity of the hypersonic waves concerned. 
Above the A-point we take c =w,, the velocity of first sound; below the A-point we 
may take with sufficient accuracy c = ug, the velocity of second sound. In order to 
account for the observed behaviour of the expansion coefficient of helium, we 
must explain two points : (i) why (du,/0V), above the A-point is negative; (ii) why 
(du,/0V), below the A-point is very small and positive. 

To explain (i) we recall the relation u?=x*/p, where y* is the adiabatic com- 
pressibility and p the density. Assuming the existence of short range forces, y* 
must depend on volume much more sensitively than p, so that the ratio will decrease 
with increase in volume. 

To account for (ii), which has recently been checked experimentally, we recall 
the relation u53=(1—p/p,)/d(1/S)/d7T. It is unlikely that either the fraction of 
normal fluid p,/p or the entropy per unit mass would be sensitive to volume 
changes, so that a very small value of (du,/0V), is to be expected. Quantitatively, 
we may investigate the magnitude of (du,/0V), by finding a relation between the 
change in velocity due to alterations in volume and the change in velocity due to 
variations in temperature. Let us write f=w,/T* =(S?p,/Cp,)?, so that f depends - 
only on the quantities which will be given by a theory of the Debye type, f=/(T/6). 
The A-point 7) is given by the condition that all atoms are normal, i.e. by 
py(T'y)=p. We have seen that p,ocH/w3. Hence Ty will be given linearly in 
terms of 6, and, therefore, f=f(7/To). Vhus (0f/0T)) = —(T/T,)(Af/eT) and 


hence (x ~ <P ydinT, af 
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The slope of the A-line in the p, T diagram of helium is negative, so that 
(dln Ty)/dIn V),» is positive, although very small. The sign of (0f/0T),, i.e. the 
slope of the curve of u,/T? against temperature, is found to be negative in the tem- 
perature region so far investigated. Hence (du,/0V),=T?(0f/0V), will be very 
small and positive, as required to explain (ii). In principle, this argument might 
be extended by the methods of §5, so as to become independent of experimental 
values of velocity. 
o7; THE ORDINARY THERMAL CONDUCTIVITY AND THE 
ATTENUATION OF SECOND SOUND 

Debye (1914) and Brillouin (1914, 1920) have suggested a theory giving the 
coefficient of thermal conductivity of solids (and liquids) in terms of the velocity 
and attenuation of the hypersonic waves of sound supposed present. ‘The theory 
is analogous to the well-known one for the thermal conductivity of a gas, i.e. that 
leading to the result K=4mCnlc, where m is the mass of a gas molecule, n the 
number of gas molecules per unit volume, C the specific heat per gramme, / the 
mean free path of a molecule and c the mean velocity of the molecules. (mCn6T) 
is the heat transported across a temperature difference 5T by the molecules con- 
tained in a unit volume; the analogous quantity in the wave picture of a solid or 
liquid will be 
cee a 
OT \ c& exp(hy/kRT)-1 
per unit frequency range, since the bracketed expression is the thermal energy per 
unit volume attemperature 7. Inthe wave picture / will now refer to the mean free 
path of a wave—.e. to the reciprocal of the attenuation coefficient of the sound 
waves concerned, and c will refer to the mean velocity of these waves. Integrating 
over all frequencies, we find for the coefficient of thermal conductivity the expression 
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For the ideal case in which c, v,, and / are independent of temperature this reduces 
to 


K=tel 
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OT), cc exp(hv/RT)-1 Y 
where C is the specific heat per unit mass. This ideal relation is the same for - 
solids, liquids and gases, provided the symbols are interpreted in the appropriate 
manner. ; 

For the case of liquid helium we have K =4p(,u,C, +/,u,C,), where C, and C, 
are the contributions to the specific heat from first and second sound respectively. 
As we have seen, below the -point u,C, > u,C,, whilst above wu, C, =0, so that we 
have 

K=~4pl,u,C(T) for 0<T<)-point, K=4pl,u,C(T) for T>A-point. 
The value of the product uC does not vary by more than an order of magnitude for 
liquid helium over the greater part of its temperature range. ‘This is due to a 
series of compensating effects. For instance, the large value of the specific heat 
just below the A-point is compensated by a fall in the velocity of second sound, and 
the change over to the value of the velocity of first sound is compensated by the 
fact that the specific heat just above the A-point is much smaller than that just below 
it. Thus any striking variation in K would be due to a change in the attenuation. 


e 
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We must stress the fact that K is the ordinary thermal conductivity, 1.e. the 
heat flux per unit temperature gradient. As has been shown (Dingle 1948), a 
finite value for K will result in absorption of second sound by irreversible processes. 
If we assume that the attenuation due to this cause is of prime importance, we may 
find J and K explicitly by combining the results of this section with a result of the 
previous paper, i.e. /, =2u3pC/Kw?, where us, p, C and K refer respectively to the 
velocity of second sound, the fluid density, the specific heat per unit mass and the 
ordinary conductivity, and w/27 is the frequency. This expression shows that the 
mean free path is in fact a function of the frequency v. Strictly speaking, therefore, 
we are unable to use the simplified expressions giving K in terms of /, except for 
rough estimates. We may show very simply that the theory gives reasonable 
value for K and 7. The minimum free path J,,;, will be that for the waves of 
highest frequency v,,. Thus 


ee fie 87uspC/Kv;, , and K ard pl intl 
giving K ~ (812/3)!pu2C |v, =(872/3)! npu2C/R8, 
and Dnin ~ (24)? 21g] Vp, = (24)? hug/ RO. 


Numerically, J,,;,-~7 x 10-*cm. at 1:5°K. This is a satisfactory value, showing 
that the minimum free path is of the order of twenty interatomic distances. A 
smaller value would have shown that the maximum frequency of vibration is 
dictated not by the number of normal modes available, but by absorption effects. 

We may proceed further and determine the approximate mean free path of 
second sound waves of attainable frequencies. Since U(v)a1/v? for a given 
temperature, we have [(v) ~7 x 10-*(v,,/v)? at 1-5°K., where v,,~4 x 10! cycles 
per second. 


S35) DHE DEBE DHE ORY VANID) Ate NE @r Si Sileinvan bh OR 
MODIFICATION 


In this section we wish to show that the use of the ordinary Debye theory is not 
justified if the sound velocities (at constant volume) vary with temperature. 

The number of independent modes of vibration per unit volume is known to be 
47rv?/c3 per unit frequency range, where c is a velocity of sound. We shall find it 
convenient to rewrite this expression. Let p be the momentum of a wave, then 
the energy is given by « =pc = hy, so that the number of modes per unit volume is 
given by 47p?/h? per unit momentum range. The upper limit to this momentum, 
Pm Say, is given by the consideration that the number of modes per unit volume is 
3N, where N is the number of atoms per unit volume. Thus 3N= (4n/h*)| a pide, 
and hence p,,=(9Nh*/47)'. Therefore p,, is independent of vemperaine 
However, the energy of each wave of given momentum, ¢ = pc, may not be constant, 
since the velocity c may vary with temperature. Supposing for the present that 
this variation of the velocity is negligible, we have for the partition function 
Z=I1;{1 —exp(pc/RT)}, where p; are the allowed momenta, this being the well 
known result valid for a set of oscillators with no restriction either on the number of 
quanta allowed per oscillator or the number allowed per energy level. 

We may proceed in two ways: either we may work out the free energy 
A=—kTI\nZ and from this the thermodynamic quantities, or else we may find 
the mean occupation numbers ; from the relation 7,;=2,01ln Z/dz,, where 
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%;=exp (—p,c/RT), and from these the mean energy, and hence, by differentiation, 
the specific heat. Using the former method, and still supposing the temperature 
variation of the velocity to be negligible, we have 

A=—-kTInZ=kT2,In {1 —exp(—p,c/RT)} 
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Introducing the Debye characteristic temperature @ by the relation R#=p,,c, we 
obtain 
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per molecular volume. Using the well-known thermodynamic formulae for the 
entropy S, the energy E, and the specific heat C, viz. S=(E—A)/T=—(@A/0T), 
and E= —T?[0(A/T)/0T], (from the Gibbs-Helmholtz equation), and 
C=—T(@A/oT’), =F /eT, 

we deduce equations (1), (2) and (3) of §3. 

Combining (8) with the relationship p= —(@A/0V),, we find the equation of 
State 

paw idV=—Eldinvsov |p. 6 0 Fens (9) 

where W is the energy at the absolute zero. 

Using the second method, that involving occupation numbers, we have 

n;= 2,0 In Z/0z;=1/{exp (p,c/RT) — 1}, 
so that «;=p,c/{exp (p,c/RT) —1}, and, therefore, — 
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which is the same as the expression deduced by the first method. It is clear, 
however, that the two methods would yield quite different results if we took into 
account the variation of the velocity with temperature, for then the value of / given 
by the second method would be unchanged, since no differentiations with respect 
to temperature are involved, whilst the energy derived by the first method would 
be different. 


§9. TEMPERATURE-DEPENDENT DISTRIBUTIONS OF 
ENERGY LEVEES 

The difficulty pointed out in § 8 is caused by the fact that the energy levels which 
we are considering are themselves functions of temperature. ‘This nay be seen 
from the energy relation «=pc, where p may be considered as independent of 
temperature, since we must integrate over all values of p up to a maximum value /,, 
which is independent of temperature. The concept of temperature-dependent 
energy levels arises only because of the difficulties of a complete treatment by the 
methods of statistical mechanics. If we included in the partition function every 
possible arrangement of the system, temperature-dependent energy levels would 
never arise. However, in many problems which are treated approximately, 
temperature-dependent energy levels are bound to be the rule. It is known in 
several special cases in classical statistical mechanics (or rather quantum statistical 
mechanics at higher temperatures) that the correct procedure is. to replace the 
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energies in the exponential terms by free energies. Rushbrooke (1940) has 
considered this problem in more detail, and has pointed out that we obtain 
agreement between statistical mechanics and thermodynamics by this procedure 
in any case in which, in the notation of his paper, N=A(0/0A)In G(Aé*) and 
E =6(0/00) In G(A6*), where Ind = —1/RT and G is a generating function: A is 
found by equating the sums of terms suchas N to the total number of systems in the 
assembly, and 6 by equating the sum of terms such as F to the total energy of the 
assembly. ‘This is equivalent to the statement that the procedure of replacing the 
energy by the free energy in the exponential terms is justified for all cases in which 
the method of steepest descents is applicable to the contour integral appearing in 
the evaluation of the partition function. That the Beg is one of great 
generality is shown by the following demonstration. 
In the customary theory, the partition function is taken to be 


Z=Xexp|[— Ln, <,/RT], 


where the «; designate the energy levels and n, the occupation numbers. ‘The 
outer summation is carried out over all possible configurations (e.g. in the case of a 
fixed number of particles, N say, this is subject to the condition that U,n,;=N), 
and the inner summation according to ell configurations consistent with the parti- 
cular statisticsemployed. Letus replace this formally by Z=Xexp[—,n,f, RT]. 
From the Gibbs-Helmholtz equation A—E=T(dA/0T),, where A is the free 
energy, we obtain H= —7?(0(A/T)/0T],. By definition, A= —kRT In Z, so that 


E= nr (2) _ Dial fi— TOF/OT)v}] exp [— Ein, fi RT] 
oT Jy Nexp[—X,f,/RT] 

2 [2 n,<;] exp[— U7fi/RT] 

Lexp[— dn, f,/RT]  ’ 


Clearly we must have F= 


since E is merely the energy of each configuration, Un,«,, averaged over all possible 
configurations. Hence, in general, 


e=fi— TefleT)y 


which is the Gibbs-Helmholtz equation for the energy levels, provided the f; are 
interpreted as free energies.. 

This result may easily be confirmed for the particular partition function in 
which we are interested, i.e. Z=I1,{1—exp(—f,/RT)}—?. 

The effective free energy of a sound wave may be found from tke expression 


<=pc(T). From the eiaheiaiee eo equation, /-- T(@f/0T),;,=.«, we have 
[oe (FDOT |p <</T? and thus f=7 [t «/T*)dT, where the upper limit 7, is left 


undefined. In the sequel we shall inl write c for [c(T)]p. 
We take for the logarithm of the partition function the expression 
Aq (Pm 


nD 
lie BB, p*dpin| 1 exp} ~ AIe celryart |, 


T, 
where p,,=h(9N/47)'. For convenience we put x= (p/R) | (c{T*) dT. so that 
in 


h (9N\ (7 ON f: 
X= alge) ie mae: henoin 7 a , ” dein (1 —e*) and thus 


m 
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A ORs fa 4 Ae ered ae 
Baa kinZ=— | © stdin (1 *)=3R{ In(1~e ere aS 
Sree (9) 
per molecular volume. ‘The energy E may be obtained either directly, thus 
cy 4a far (pc)p* dp a 9Re il ie ax, 
he )o expl(p/k) | (clT)4T]-1 [7° (eT) aT *m io OT 
Bo in ies halal RAG acca (10) 


or from A/T by differentiation, noting that 
1 OX, a4 c/T? 


ToD if (c/T2)dT 


‘The solutions given by the two methods now agree. 


The entropy may be found either by using the thermodynamic definition 
A=E—TS, giving 


— tm x3 
gs ~-4 _3p In(1-e™) + = | mee it ee emle aay 
1 Soo eek Tee (c/T?) aT | | 


Raters (11) 

-or from the thermodynamic formula S= —(0A/0T),. There is no need to verify 
that we obtain the same result, since by elimination of S between these two 
expressions, i.e. S=(H—A)/T= —(0A/0T),, we return to the Gibbs-Helmholtz 
equation which was used in the derivation of the energy : the two values of S are 
therefore bound to be identical. Finally, the specific heat C may be found either 
by differentiation of the energy, or of the entropy, i.e. C=0E/0T=T(OS/0T),. By 
-either method we obtain 


hy (Cun anx c 2 GG ans 
ele oak Jlesra * Reimar | 
K \ L / ; J 


Ly 
T 


Dales oP Aa Ge 7 as pee Nee a ei 


In a similar way, we may show that the equation of state is now 


dis ey Ma ke 
aM 


P+ a a cas = ee (13) 
-and that the equation p,/p should now read 
4 E 
Pap = ea A i (14) 
EG ear: 


The theory given above leaves undetermined the upper limit 7) of the integrals 
involved. In our case three choices might be made: (i) 7)=0°K.; in this case 
we must assume that c(7) varies in such a way that the integrals do not 
‘diverge at extremely low temperatures, i.e. near T=0, c(T) «T”, where n>1; 
Hat) 7, =2:19°xK,, the A-point; (11) T5= ©. 

The second and third give the same result, since u,=0 for T>2-19°K. 
However, the integrals would have a zero value at the A-point, which would lead 
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to difficulties—for instance, infinities would appear unless the attenuation of the 
sound were taken into account in detail. It is therefore probable that T) =0°K.., 
predicting a zero value for the velocity of second sound at the absolute zero. 


Note added in proof. Recent Experimental Work. 


Peshkov (1948 a) has recently made measurements of the velocity of second 
sound at lower temperatures than those used in previous experiments. These 
show that the velocity again rises at temperatures lower than 1:12°K., in 
accordance with Landau’s prediction (Landau 1941, Dingle 1948, §11). 
Thus the mean value of 0, will be greater than that estimated in §5 of the 
present paper. As a result, the values of the thermodynamical quantities will 
be reduced, leading to better agreement between theory and experiment. 

Peshkov (1948b) has also measured the velocity of second sound as a 
function of pressure, i.e. (@u/0p)p. This he finds to be small and negative, 
so that (@u,/0V), will be small and positive, as predicted theoretically in § 6 of 
this paper. Quantitatively, by similar arguments to those of §6, we obtain 


Ouy\ — TdT, (Oug/V/T) 
(3) ECP), 

This relation yields values of (0u,/0p), in close agreement with Peshkov’s. 
experimental results for temperatures near the A-point—. e. in the temperature 
region in which wy varies rapidly with temperature. For temperatures in the 
neighbourhood of the maximum of u,(7) it gives values which are too small by 
a factor of about four. It follows from the above relation that —(du./0p)>- 
should again rise in value for temperatures lower than 1-12°x. 
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The Surface Energies of Inert-gas and Ionic Crystals 
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ABSTRACT. It is shown that the surface energies of inert-gas crystals are proportional. 
to their sublimation energies. "The constant of proportionality is evaluated for the {111},. 
{100} and {110} faces for different assumed forms of the interatomic repulsive forces; to 
within 15°% it is independent of the face, and of the form of the repulsive force. The 
distance of the outermost plane of atoms from the next at a {100} face is 2:54°% greater 
than normal; this distortion makes a negative contribution of 1% to the surface energy. 

The surface energies of sodium and potassium halide crystals are calculated; the ratio: 
of the energies of {110} and {100} faces is about 2:5. The contribution of van der Waals. 
type forces to the surface energy of a {100} face is almost half the total. 

Elementary methods are given for the evaluation of lattice sums. 


§4-INTRODUCTION 


HE total free energy of a volume of liquid is the sum of two terms, one 
proportional to the volume, the other to the surface area; the surface free: 
energy per unit area is defined as the proportionality factor in the second 

term. ‘The surface free energy per unit area of a crystal fa-e is similarly defined, 
having regard to the fact that different faces of a crystal have different energies. 

Both the surface free energy and the surface “‘internal”’ energy are difficult to. 
calculate because they depend on the thermal vibrations. ‘The surface energy at 
0° k. is asimpler quantity, for except for zero-point vibrations the molecules are at: 
rest. Ifthe molecules can be regarded as acting upon each other with forces that. 
depend only upon their separation, then the surface energy at 0°x. is the excess. 
potential energy of molecules near the surface. 

An increase in area of a given volume of material can occur by its division into 
two parts by a plane, and the separation of these parts to an infinite distance. ‘This. 
can be supposed to occur in two stages and the energy change of each stage can be: 
calculated. In the first stage the molecules of each part are supposed held in their 
original relative positions by constraints, so that the energy per unit area needed to. 
separate these two parts to infinity is 2Uj, where —2U) was the mutual potential. 
per unit area between the two parts before separation. In the second stage the 
constraints on the molecules are relaxed, the surface molecules are no longer in 
equilibrium and move to a position of lower energy; this reduces the energy of 
each surface by an amount Uj per unit area. Thus the surface energy per unit. 


area at 0°K. is 
OE OA OF a re eee eA (1) 


Uj is a small correction to U5. 
Numerous authors have attempted to calculate Uj for homopolar liquids on 
the assumptions that the forces between any two molecules act along their line of 
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centres, are independent of the presence of other molecules, and depend only on 
their distance apart, i.e. they are central and additive. Using these assumptions 
Gauss (see for example Rayleigh (1890) or Fowler and Guggenheim (1939)) 
treated a liquid as a continuum and obtained the expression 


oh?) 
y= — den? | WEN GR | Aq « eee (2) 

YO 
where e(r) is the potential between two molecules at a distance r apart and m is the 
number of molecules per unit volume. This gives correctly the contribution to 
U,’ of those molecules distant from the plane of section.* » It shows that if €(7) 
decreases less rapidly than r~ the integral does not converge, and then the concept 
of a constant surface energy per unit area independent of size is not valid. Ifa 
volume is so large that gravitational forces between the molecules affect capillary 
phenomena, they must be estimated independently. 

If, however, <(7) contains terms that, at y=0, change more rapidly than r~4, 
the integral (2) diverges at the lower limit. Edser (1922), Bakker (1928), 
Bradley (1931), Kasseland Muscat (1932) and Belton and Evans (1941) have avoided 
this divergence by using the atomic diameter as the lower limit; this is quite 
arbitrary and introduces very large and unpredictable errors. 

To obtain an accurate estimate of Uj it is not sufficient to treat a liquid as a 
continuum; a knowledge of the molecular arrangement is necessary. Jura (1948) 
used the experimentally determined radial distribution function of argon and 
calculated Uj, p, the surface potential energy of the liquid at the melting point, 
from a generalization of equation (2). He assumed that the distribution of atoms 
near the surface is not altered by the presence of the surface. ‘This procedure is 
somewhat tedious and demands a knowledge of the radial distribution function. 

In this paper it is shown that for a homopolar substance, either liquid or solid, 
in which the interatomic forces are central and additive, 


U, Kee, | 3s a eee (3) 


where ey is the heat of evaporation and ay is a characteristic interatomic distance; 
both quantities are supposed extrapolated to 0°K. 'The value of K depends on 
the form of the interatomic force, <(r), and upon the molecular arrangement. 
Values of K for {111}, {100}, {110} faces of face-centred-cubic crystals and for 
different interatomic repulsive forces are calculated by a method of lattice sums 
devised by Born and Stern (1919) for ionic crystals. The theory is applied to 
inert-gas atoms, for which the assumptions are most likely to be valid. 

Calculations of the surface energies of the {100} and {110} faces of sodium and 
potassium halides have been made by Born and Stern (1919), Lennard-Jones and 
‘Taylor (1925), and Dent (1929). However, doubt has been thrown on these 
results because they were less than half the observed value for the liquid. It is 
shown that this discrepancy is largely due to the neglect of van der Waals attractive 
forces, which in these ionic salts make a large contribution to the surface energy. 

No reliable measurements have been made of the surface energies of solids and 
it is necessary to compare the theoretical results with the observed values for the 
liquids. _ Exact agreement is not to be expected because of the different molecular 
arrangement in solids and liquids, especially at the surface. 


_ *'The contributions of distant molecules to Uj are found by multiplying 7 of equation (30) 
‘by—4hn; for a continuum h->0 and the second integral vanishes. 
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§2. THE MUTUAL POTENTIAL BETWEEN TWO PARTS OF 
AIG RAYS) Age, 

It is necessary to calculate the mutual potential per unit area between two parts 
of a homopolar crystal separated by an imaginary plane of section. Because of the 
regular crystal structure all atoms can be regarded as lying on planes parallel to the 
plane of section, and when all lattice points are equivalent the separation h between. 
all adjacent planes is the same. The mutual potential between the two parts, 
Land II of the crystals (see Figure), is the sum of the potentials between all planes in 


ee ey ER 


0 
Swern S oe ee Ol 
ee ee 


part I and all planesin part IJ. The potential between any two such parallel planes 
depends only upon their separation. 

From the Figure it is seen that in the expression for the mutual potential the 
contribution of a pair of planes at separation h occurs once, at separation 2h twice, 


and so on. If the planes in part I be enumerated /=1, 2, 3,...., then the mutual 
potential (per atom in the surface plane) —2u, is 
Bey SER fer ygun = athe) tla 204i had (4) 


i>1 
Here uj is the surface energy per atom in the surface plane, r is the distance 
between an atom O in the surface plane (/=0) of part Il and an atom in the plane / of 
part I; the summation is over all atoms in part I. 

This method of calculating the surface energy is based on the assumptions that 
the interatomic forces are central and additive; these are valid only for inert-gas. 
and ionic crystals. In inert-gas crystals the attractive forces between atoms are 
van der Waals, which London (1937) showed to be additive. The attractive 
forces between the ions of ionic crystals are van der Waals and electrostatic forces ; 
the latter are clearly both central and additive. For neither inert-gas nor ionic 
crystals is it certain that the repulsive forces are additive, but this will be assumed 
for the inert-gas crystals; for ionic crystals only repulsive forces between nearest 
and next-nearest neighbour ions will be considered. The form of this assumption 
is not important because it affects only the interaction between atoms which are not 
nearest neighbours, and this is small because the forces decrease rapidly with 
distance. 

The method of calculation is not appropriate to metals for which the surface 
energy is best regarded as arising from the change in kinetic energy of the free 
electrons due to the presence of the surface (Huang and Wyllie 1949). = Nor can 
it be applied immediately to molecules with a permanent dipole moment, for these 
will be oriented at the surface of the liquid, and this will give rise to large values of 
U’. For covalent crystals, such as diamond, the best model is perhaps that of 
Harkins (1942) who considered only forces between nearest-neighbour atoms ; 
this leads to a simple expression for K. 


170 R. Shuttleworth 


§3. THE SURFACE ENERGIES OF INERT-GAS CRYSTALS 
(i) Calculation of U% 


The force between two atoms is attractive at large distances and repulsive at 
‘small distances ; at intermediate distances there is an equilibrium separation where ~ 
the potential energy isaminimum. A simple expression that shows this behaviour 
ds 

(7) =Ar4 —pr st, be 2 SOR 


where «(r) is the potential between two atoms a distance 7 apart; —sAr~S* can 
be regarded as the short range repulsive force and fur~“*) the longer range 
attractive force. From a knowledge of ¢(r) both uj and uj can be calculated. 

The lattice of inert-gas crystals is face-centred cubic; this may be regarded as a 
simple cubic lattice, of lattice constant a, in which only alternate sites are occupied. 
"Taking coordinate axes parallel to the cube axes, the coordinates of an atom, with 
respect to another, are (121d, M9, M34)) where my, M9, mz are integers (positive, 
negative or zero) such that m,+m,+msiseven. Combining equations (4) and (5) 


Sie os: arn | 


11 (m* + m3 + m>)s? 
M+ Ms+ M3 eVeN 
l 
al 
— pa >» ILO. Se EOn SLO yaa 
ae | I>1 (my tty + ma) I (6) 
M+ M.+ Ms even 
/ pias, Mu" -s a =f 7 5 
—2u/,= BiAajh — Bias". 2) se wee (7) 


Numerical values of the lattice sums B” and details of their computation are given 
in the Appendix for {111}, {110}, and {110} planes for powers between 6 and 12. 

Instead of using explicit values of the force constants A and p it is possible to 
‘express My in terms of the heat of sublimation of the crystal. The mutual potential 
-of an atom in the centre of a crystal with respect to all the other atoms is 


— 2) = Xe(r), Jt eee) 
where ¢€, is the heat of sublimation per atom at 0°K. (corrected for zero point 
vibrations) and the summation is over all the other atoms. The factor 2 occurs 


because when the sum is repeated for every atom as origin the interaction of each 
pair is counted twice. Using equation (5) 


a ] 
— Leg = ay? | ~ PB Ha Meee 
M,+ms+msz even (m1; + M3 + mi) 5 


1 
= —t Ss C 
‘aad | ~ TREO LESUIN Ilo. CR ooc (9) 
M+ M+ Ms even (mi; =F ms + m3)? 
dey = AG = Alay Se eee (10) 


Lennard-Jones and Ingham (1925) have evaluated A” for powers between 4 and 30. 
At 0°K. the lattice potential energy eg is a minimum with respect to variations in a. 


cA ay = tA sags GS eee alt) 
Eliminating Aap* and jap’ from equations (7), (10), (11), 
sBy/ Al ~ tBt) A! 


— h(t | 1a) eeeee (12) 
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The attractive forces between inert-gas atoms are van der Waals type, for which 
t=6. In Table 1,A(s,6) for a {111} plane is tabulated for values of s between 7 
and infinity. Increasing values of s correspond to increasing “hardness” of 
the atoms ; when sis infinite the atomic model is that of rigid spheres with attractive 
forces. 


Table 1. Variation of ,,,2 with ‘“‘ Hardness” of the Atoms 


s 7 8 9 10 11 12 oe) 


uik(s,6) 0°502 0-450 0-418 0-398 0-384 0-374 0-313 


If the interaction between atoms that are not nearest neighbours is neglected, 
then B//A/ and Bf/A‘ are each equal to the ratio of the number of nearest neigh- 
bours in the plane /=1 (3), to the total number of nearest neighbours, (12); & is 
then +. 

The most significant basis of comparison for the surface energies of different 
crystal faces is the energy per unit area Uj 


pUp=pti/wpaj=pK(s, tea, a nee (13) 


where 1/w aj is the number of atoms per unit area in the surface plane of crystal 
face {P}. In Table 2 values of »K(s, 6) are tabulated for {111}, {100}, and {110} 
planes for s between 7 and infinity. 


Table 2. Variation. of K with Crystal Face and ‘‘Hardness”’ of the Atoms 
s 7. 8 9 10 ital We co 


nK(s,6) 0-290 0-260 0-242 0-230 0-222 0-216 0-181 
sooK(s,6) 0-292 0-264 0-247 0-236 0-229 0-224 0-196 
noK(s, 6) 0-301 0-273 0-257 0-247 0-240 0-235 0-206 


The variation of surface energy with crystal face arises from the different 
arrangement of the atoms that are near the surface; the contribution to pUj of 
atoms distant from the plane of section is the same for all faces. ‘The greater the 
relative importance of the interaction between neighbouring atoms the greater is 
the variation of surface energy with crystal face; thus the variation increases with 
increase ofsand¢. Yamada (1923, 1924) and Frenkel (1945) have shown that the 
surface energy per unit area of a crystal changes continuously with change of orienta- 
tion of a crystal face, and that faces with large indices need not necessarily have 
large energies. ‘lhisis because faces with large indices can be regarded as a series of 
steps formed from planes of low indices and of height equal to the interplanar 
separation. When such steps are few in number the increase in surface energy is 
only slight. Thus for a homopolar face-centred-cubic crystal the {110} and {111} 
faces have the absolute maximum and minimum surface energies per unit area; 
for an inert-gas crystal they do not differ by more than 10%. 

From a knowledge of the heat of sublimation and the lattice constant of an 
inert-gas crystal it is possible to calculate the surface energy to within 15%; a 
more accurate estimate demands a knowledge of s. Values of s are derived from 
the variation of the second virial coefficient of the gas with temperature; Bucking- 
ham (1938) has shown that for inert gases a form of e(r) in which t= 6, s=12 gives 
good agreement with experiment. Values of Uj for the {111}, {100} and {110} 
faces of the inert-gas crystals have been calculated from their heat of sublimation, 
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using K(6, 12), and are given in Table 3. They are compared with the value of the: 


surface energies of the liquids extrapolated to 0°K. (Guggenheim 1945). 


Table 3. Surface Energies of Inert-gas Crystals 


€o ag* 115 100U 1109 Uo iia. 
(erg/atom) (A) (erg/cm?) — (erg/cm?) (erg/cm?) -(erg/cm?): 
Ne 4-08 x 10714 2-26 Ze? 17-9 18-7 15-1 
A 13-89 2:70 41-1 42-7 44-7 36°3 
Kr 19-23 2°79 53°3 5°35 57°9 
Xe 26°87 3-09 60-7 63-0 66-0 


* The lattice constant is Ady. 


The heats of sublimation at 0°. were calculated by numerical integration: 
from the latent heats at the boiling point, the heats of fusion and the specific heats. 


of the solids and liquids ; a correction 9R@/8 for the zero point vibrations was made, 


© being the Debye temperature. Nocorrection was made for the departure of the 


gases at atmospheric pressure from the perfect gas law. The experimental data 
are those of Clusius and his collaborators (1929, 1936, 1937, 1938, 1939). 

The surface energies calculated for the crystal are greater than those found by 
experiment for the liquid. For inert gas crystals dU %/da is negative, but the 
difference in density of liquid and solid is not quite sufficient to account for the 
discrepancy. 

The effect of zero-point vibrations on the surface energy has been neglected ;. 


for while they would invalidate any calculation for helium their effect on the surface ~ 


energies of crystals formed from the heavier atoms will be small. In (ii) below it is 
shown that the decrease of energy Uj due to surface distortion is less than the 
uncertainties in the calculation of U>. 


(ii) Contribution of Surface Distortion to Surface Energy 


In inert-gas crystals the separation of the atomic planes parallel to the surface 
is greater near the surface than in the interior of the crystal. ‘The surface energy 
of the distorted crystal is less than that of the unstable undistorted crystal by U/. 

The distortion arises because of the interaction of atoms in non-adjacent 
planes. In the interior of a crystal (see Figure) the net force between all the atoms 
above the atomic plane /=1 on all the atoms below the plane is attractive. This is 
because, in homopolar crystals, the attractive forces decrease more slowly with 
distance than do the repulsive forces, and the separation of the parts of the crystal 
above and below /=1 is greater than the equilibrium distance. This attractive 
force is balanced by net repulsive forces between the atoms in the plane /=1 and 
the two parts of the crystal above and below /=1. If the part below /=1 is 
removed, so that /=1 becomes a surface plane, there is an unbalanced repulsive 
force which causes the atoms in the plane /=1 to move away from the rest of the 
crystal until the attractive and repulsive forces again become equal. Similar 
distortions occur for atoms in successive planes but these are very much smaller 
and will be neglected. ‘The distortion will be least on crystal faces formed from 
close-packed planes, for these planes are most widely spaced and the forces between 
non-adjacent planes smallest. If it is possible to neglect the interaction between 
atoms that are not nearest neighbours the surface distortion is zero. 


“ 
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U; is the decrease of potential energy when the surface plane of atoms moves 
from the position it occupied when at the centre of the crystal to its equilibrium 
position of minimum energy. Calculations of distortion and of U,” were made for 
the {100} face of an inert-gas crystal (s=12,t=6). Distortion occurs by an equal 
displacement, in a direction perpendicular to the surface, of each atom in the 
plane/=1. The potential energy per atom of the surface plane with respect to 
the remainder of the crystal is: 


p=C Ado — Couas® wacertabte Rates (14) 
where (OAC ee Oa) 7)’ cutee 2.0 C8 Oy eee es (15) 
I>2 


and cay is the distance of the surface plane, /=1, from the next plane, /=2. In 
equation (15) 7 is the distance between a particular atom in the surface plane and 
any atom in the remainder of the crystal; summation is over all atoms in the re- 
mainder of the crystal. In Table 9 values of CG and Cj, for a number of values of c 
are given. ‘lhe equilibrium position of the surface plane occurs at the value of ¢ 
for which ¢ is a minimum, viz. 


GTA DAT SAA CUdee St Aes (16) 


Aap? and pap ® have been eliminated by equation (11). A solution to equation (16) 
was found by numerical differentiation: ¢,,y;,=1-0254. Thus in an inert-gas 
crystal the distance between the surface atomic plane and the next is, for a {100} 
face, 2-5°/, greater than the distance in the interior of the crystal. 

Values of the lattice sums at the equilibrium separation ¢,,,,;, were obtained by 
interpolation and are included in Table 9. Then 


Pe ep NC ola SACWa.. | 1. ast? (17) 
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where u/, is the decrease of ¢ as c increases from 1 to 1-0254, and AC}, and AC{ are 
the corresponding decreases in the lattice sums. As in the derivation of equation 
(12) 

Mpg C ioe AAG a Aglegy 9 wares (18) 


see I SE VOY re I a oa (19) 


10U Gis less than 1°% of y99U4 and can be neglected because it is smaller than the 
inevitable errors in the calculation. Surface distortion will reduce the energy of 
a {111} face in a somewhat less, and of a {110} face in a somewhat greater propor- 
tion, but not sufficiently to affect appreciably the relative magnitudes of different 
faces. 


M4 DHE oUREACE ENERGIES OF TONIC CRYSTALS 
(1) Calculation of U% 


In crystals of the alkali-halides the alkali and halide ions are arranged on 
alternate points of a simple-cubic lattice. The halide ions have a negative charge 
of one electron and the alkali ions an equal positive charge. ‘There are also van der 
Waals attractive forces and short range repulsive forces between the ions. Mayer 
and his colleagues (Born and Mayer 1932, Mayer and Helmholz 1932, Mayer 1933, 
Huggins and Mayer 1933) have considered these forces in detail and given numeri- 
cal values for the force constants. According to them the potential between two 
alkali ions at a distance r apart can be written 


ei (r)=e% tc, 78d 7? +1-25b) exp(—r/p);  ...... (20) 
PROC. PHYS. SOC. LXII, 3-——A 13 
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the potential between two halide ions 

e__(r) =e"! —c_-r-*—d__r-*+0-75b2-exp(—r/p), —...... (21) 
and the potential between an alkali and halide ion 

é4 (n= — 67 =¢, 7 dF 5-0 eXD| 9) B)ae ae oie (22) 


The van der Waals terms proportional to the inverse sixth and inverse eighth 
power of the distance are referred to as the dipole-dipole and dipole-quadrupole 
interactions respectively; the latter term is small compared to the first. ‘The 
repulsive potentials decrease exponentially with distance ; in a crystal it is necessary 
to consider only the repulsion between an ion and its six nearest neighbours, 
which have opposite charge, and twelve next-nearest neighbours, which have the 
same charge. 

If an ion is selected as origin all ions with the same charge have coordinates for 
which mm, +m +mz is even, ions with opposite charge have coordinates for which 
mM, +my+ntz is odd. All planes (except {111} planes) contain equal numbers of 
alkali and halide ions. To find the mutual potential between two parts of an 
alkali-halide crystal it is necessary to consider the potentials between (a) the 
alkali ions of part I and the alkali ions of part II, (0) the halide ions of part I and the 
halide ions of part II, (c) the alkali ions of part I and the halide ions of part II, 
(d) the halide ions of part I and the alkali ions of part IJ. The positions of the 
alkali and halide ions in the lattice are equivalent so that (c) and (d) are equal. 
The electrostatic contributions if considered separately diverge, and they must be 
grouped in pairs (a) and (c), (6) and (d), to obtain convergence; these two pairs 
are equal. As in equation (7) the mutual potential per alkali (or halide) ion in the 
surface plane is 


—2u) = —2B,' e? ap! ~ Bole. +¢__)ap® —2Boc,_ah® 
— BY(d,.,+d__)ay§ —2Bid,_ay® 
+ .N"(1:2562. +.0:75b7) exp (— 1/2aq/p) 


An ZNO b exp (— Gale), ay See ee ee (23) 
where dy is the alkali-halide ion distance, B/”’ is an electrostatic lattice sum 
ms M+ M2+ Ms 
Bem Dy, ee I cy: | (24) 


i>1 (mi? + m3 +m)’ 

M+ M.+ ms, odd and even 
the factor (—1)'"*""*" arising because of the opposite sign of ions at alternate 
lattice points ; the B;s represent the lattice sums for the van der Waals interaction 
between alkali and halide ions: 


B= ee ee 
1>0 (my + mz + m3)" 
M,+m.+ ms Odd 


pre: (25) 


N' and N" are respectively the number of nearest and next-nearest bonds per 
surface ion that cross the plane of section. Values of the lattice sums (divided by 
2wp) for {100} and {110} faces are given in Table 4, and their method of calculation 
is described in §5. 

The surface energy per unit area, U4, is up/wpa>, where 1/wpa? is the number 
of alkali (or halide) ions per unit area in the surface plane. The force constants 
given by Huggins and Mayer (1933) for sodium and potassium halides were used 
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to calculate the contributions to Uj of the electrostatic, dipole-dipole, dipole— 
quadrupole and repulsive forces. These contributions and their total are given in 
‘Table 5 for a {100} face. 

In Table 6 are given the surface energies ,,)Uj of a {110} face, the ratio 
1106/1909 U5 and a). For consistency a) was calculated from the force constants. 

The electrostatic term is predominant in the lattice energy, and the repulsive 
and van der Waals terms are only one-tenth and one-thirtieth of the electrostatic 
term. ‘This is not so for the surface energy : here the electrostatic term (positive) 
is only a little greater than the repulsive term (negative), so that van der Waals 
forces give rise to almost half of the net surface energy. This is because U{ is the 


‘Table 4. Lattice Sums used to Calculate the Surface Energies of Alkali-halides 
BY /2ap Bi/2op B,/2op Bi/2wp B,/2op N"/2wp N’'/2ap 


{100} 0-01631 0:17670 0:33820 0-06920 0-26898 1:00000 = 0-25000 
{110} 0:03174 0:-18618 0-43482 0-07353 0-36972 1-06066 0-35355 


Table 5. The Contributions of the Different Forces to U{ for 
{100} Faces (erg/cm?) 
NaF NaCl NaBr Nal KF KC] KBr KI 


Hlectro- \ 618-0 345-6 289-7 227-8 396-5 -- 250-3. 215-8 175-8 
static : 


Dipole— \ 82:8 79-6 77-9 76-2 84-1 69-8 65-0 60-6 
dipole 

Dipole— 
sa 8-8 8-6 8-4 8-4 8-3 fs 6-6 6:3 
pole 

Repulsives — 500-3 —245-3- —199°3- 151-9 294-3 — 164-3" —136-2 106-4 
Total 209 189 177 161 195 163 151 136 


Table 6. The Surface Energy, U(, of a {110} Face 


NaF NaCl NaBr Nal KF Kel) KBr ‘KI 
a10U, (erg/cm?) 640°5 445-3 396-2 SN) EAA Rea | SSeS SO, 273-6 
1100 o/100U 4 3-06 2°36 2-24. aA 2-48 2:16 209 2:01 
Qo (A.) DENY Pe — PHASES) SOA eae 05) 5095 Say? 3-482 


energy necessary to divide a crystal into two neutral parts, whilst the lattice energy 
is the energy necessary to divide it into oppositely chargedions. ‘The large ratio of 
the surface energy of a {110} to a {100} face arises mainly from the difference 
between the electrostatic terms; this is clear from Table 4. 

In the alkali-halide crystals the {111} planes are formed alternately from alkali 
and halide ions, and there appears no simple method of evaluating the electrostatic 
contribution to the surface energy of a {111} face. 


(ii) Contribution of Surface Distortion to Surface Energy 


Lennard-Jones and Dent (1928a) and Dent (1929) showed that in ionic 
crystals the most important cause of surface distortion and corresponding decrease 
of surface energy U/ is the polarization of the surface ions in the electrostatic field. 


nye 
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Verwey (1946) has made detailed calculations and shows that, in general, the 
positive ions at the surface are displaced away from the crystal and the negative 
surface ions displaced towards the crystal. He shows that for NaBr U¢ is 
32erg/cm?. It is supposed that NaBr is typical, and 19)U for the other crystals 
obtained by reducing j9)U in the same proportion. In ‘Table 7 values of i) U¢ 
are compared with the experimental values of the liquid, which have been extra- 


polated to 0°. by Jaeger (1917). 


Table 7. Surface Energies j9)U, of the {100} Face (erg/cm?) 
NaF NaCl NaBr Nal KF KCis (KBs KI 


mae i714" “455 F458 V3 ice ie ee 
Uoua) 335. 1190. «A 7Bie BSS 4s ete 


Although for alkali-halide crystals dU(/da is positive, the decrease of density 
on fusion is not sufficient to account for the discrepancy between the two values. 


Son DHE EVA LAT T ONTO SAT LICE ss UNS 
(1) Computation of B; and B; 
In order to evaluate lattice sums of the form 


l 
Bis Dy pe (26) 
t 2 See he | LD reecects 
I>1 (m2 + m2 + m3)l2 
My+M2+ Ms EVEN 
it is necessary to express /, the ordinal of the plane in which an atom lies, in terms of 


the coordinates of the atom. ‘The relation is different for different planes 
bap yk iG tg) 2, ge =U aan St ee (27) 


where m1, Mg, mz are integers (positive, negative or zero). ‘The values of / for any 
value of (r/a 9)? =m, +m +m are readily enumerated, and the contributions to 
Bj of all atoms (of part I of the crystal) that are within or on a hemisphere centre O 
and radius 7, can be found by direct summation. When f¢ is greater than 10 the 
sum converges rapidly and.an accurate value of B/ can be obtained if the contri- 
butions of atoms outside a hemisphere radius \/30a) be neglected. For smaller 
values of t the contributions even of those atoms that are outside a hemisphere 
radius 1/50a, cannot be neglected, and they must be estimated by integration. 
Before finding the contributions to B/ of those atoms outside a hemisphere 
radius 7,, a more general expression that gives the contributions of these atoms to. 
the mutual potential, —2U, in terms of ¢(r) will be derived. For distant atoms. 
each atomic plane can be regarded as a uniform lamella which extends a distance 
h/2 above and below the plane; their contribution to —2Uj is, by equation (4), 


ee [H zelda, * in ee (28) 


nis the number of atoms per unit volume, the element of volume dV is at a distance: 
r from O, and x from the plane /=0 (see Figure). The continuum representing 
part I of the crystal extends to the plane /=4; so the region of integration is that 
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portion of the semi-infinite volume bounded by the plane /= $ that is outside the 
hemisphere radius r;._ This region is the difference of two other regions: (a) the 
semi-infinite volume above the plane /=0 and outside the hemisphere rj, (b) the 
portion of the infinite slab lying between the planes / = 0 and /= 4 that is outside the 
cylinder radius 7;._ The integrals over the two regions are written in spherical 
and cylindrical coordinates respectively, and in the second integral the approxi- 
mation is made of putting the cylindrical radial coordinate equal to r, the distance 
of a volume element from O. 


ne {2 p2Qr : n (2 chia 2a 
jee an | - re(r) cos 0 sin 0 dd dO dr — oe ie \. re(r)x db dz dr. 
Bh ces (29) 
fh [ 7 fice 
I=n5 | ; e(r) dr — 7 nh | f rear. wets (30) 
In a face-centred-cubic lattice 
n= 1/205 and MiNi Gipd., oi) TE Pats (31) 


so if <(7) contains a term Ar the contribution of the integrals to B/ is 


TWp 7 


A(t—4)(ri/ao)4 — 4arp(t—2) (14/0)? 


In order that the contributions to B/ of atoms on the surface of the hemisphere 7, 
should not be included in both the sum and the integral, the lower limit of integra- 
tion 7, was defined by 


RO Need) vir 20 lhe Qa Ae (32) 


where JN is the number of atoms that are in and ona sphere centre O and radius 7,, 
and 47rj?/3 is the volume of continuum these N atoms would occupy. To 
compute Bj it was necessary to find by direct summation the contributions of all 
atoms within a hemisphere radius \/50a).___B/' is tabulated in Table 8 for {111}, 
{100}, and {110} faces for ¢ between 6 and 12. 


Table 8. The Lattice Sum B; 


t 6 7 8 9 10 11 2 
renee 056593 0:332604  0-21428 0-14391 0:09873 0:068545 0:047914 
welt. 0:70680 0:-42387 0:27680 0-18765 0-12960 0:090387 0:063389 
a Sy 1:0532 063547 0:-41596 0:28218 0-:19486 0:135863 0:-095244 


Similar computations were made for B; and the results, in the form B//2w», 
are given in Table 4. 


(11) Computation of Cg and Ci, 


In the computation of C¢ the contributions of all atoms inside and on a cylinder 
whose radius was 1/50a), and whose bases were on the planes /= 2-7, were found 
by direct summation. The contributions of those atoms in the planes /=2—7 
that are outside this cylinder were found by integrating separately for each | 
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plane; a lower limit of integration given by an equation analogous to (32) was 
used. The contributions of atoms in the planes />7 were found by a single 
integration; the lower limit of integration was taken as /=7-5. Cj, converged so 
rapidly that it was sufficient merely to sum atoms inside and on a cylinder radius 
4/14a, that lay between the planes /=2 and /=4. Cg and Cy are given in 
Table 9 for values of c between 1-00 and 1-04. 


Table9. The Lattice Sums Cg and Cj, * 


c 1:00 1-01 1-02 1-03 1-04 1-0254 
CE 0-612195 0:595823 0579923 0564479 0549472 0-571527 
Cie 0:063028 0-059367 0-055920 0-052673 0-049616 0-054144 


* JTndividual values are not accurate to this number of figures but the set is consistent. 


(iii) Computation of By’ 


Previous calculations of B/’’ for {100} and {110} faces of soditum-chloride type 
crystals have been made by Born and Stern (1919). The atomic planes parallel to 
these faces are neutral and their potentials decrease very rapidly in a normal 
direction. Using Madelung’s method (see for example Lennard-Jones and Dent 
1928 b), Born and Stern evaluated the potentials of the planes /=1, 2, 3, and found 
they could neglect the others. Their value of ,,)B(”’ is used in this paper. 

It is possible rapidly to evaluate 19)Bi/’ for a {100} plane by an elementary 
method of direct summation similar to one suggested by Slater (1939) for the 
evaluation of Madelung’s constant. If the contributions of positive and negative 
ions to B/’’ are summed separately, both diverge ; convergence is possible only if the 
net charge on all the ions included in the sum is zero. If rapid convergence is to 
occur the domain of summation must be such that all ions outside it can be grouped 
so that every group has zero charge and makes only a small contribution to the 
potential at O (see Figure). 

First, an approximate value of the potential at O due to each of the planes / = 1—5 
was found by summing over all the ions of a plane inside and on a square with centre 
on the normal cube axis through O and sides parallel to the other cube axes. In 
order that such a square be neutral an ion on a side was treated as only half in the 
square and an ion at the corner as only a quarter in the square. ‘The potentials of 
each plane /,e/do, /,@/a, 1,e/ay etc. thus calculated, converged somewhat slowly as 
the size of the square was increased. ‘The errors, which arise because of ions not 
included in the sum, have opposite signs for adjacent planes, since in adjacent 
planes positive and negative ions interchange their positions. If the planes are 
grouped in threes and the middle plane given double weight, the potential of the 
triplet rapidly converges; this is because the error of the centre sum is almost the 
mean of the errors of the exterior sums and has opposite sign. 

The two triplets, /, +2/,+/, and 2(/; +2/,+1;), were evaluated using a square 
of side 10a, ; this gave a value of Bi” (=1,+21,+3/,....) equal to 0-06525. 
The second triplet was only 1°, of the first. 
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On the Surface Free Energy of Certain Metals 
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ABSTRACT. For the purpose of calculating the surface energy and surface double 
layer for a number of monovalent metals, a simple and yet adequate model is proposed, 
which does not have the disadvantages of models employed by previous authors. In this 
the valence electrons are treated according to a simple Sommerfeld model, and the positive 
charge of the ions is distributed uniformly through the metal. It is shown that by a suitable 
combination of the methods of Thomas—Fermi and of direct wave mechanics, both the 
surface energy and the strength of the double layer can be calculated. The surface double 
layers for the metals considered are found to lie in the range 0-1-0-5 ev. The contribution 
to the surface energy from the electrostatic energy in the double layer is found to be very 
small. In estimating the surface energy, the electrons in a metal can be considered as 
moving in a box with finite potential walls of suitably determined height ¢. A temperature- 
dependent term in the surface free energy due to the change in lattice vibrations caused 
by the free surface is also considered. The surface energies obtained from ¢ determined 
in two independent ways agree closely and are found to agree with the experimental values 
of the surface tension of the liquid metals to within a factor 1°5. 


$45 INTRODUCE ION 
N spite of the high values of the surface tension exhibited by metals, it is an 
interesting fact that the ratio of the surface energy per surface atom to the 

heat of vaporization per atom is for monovalent metals of the order of 1 :7, 
for zinc and cadmium about 1 :3-5 and for mercury about 1 :2-4, whereas for 
many other substances the ratio is of the order 1 : 2. 

Shuttleworth (1949) has discussed the values of this ratio arising from various 
simple laws of force, applicable in ionic and molecular crystals, and finds results 
varying from 1 :2 to 1 :4 according to the law of force. 

In the present work we show that the peculiarly low values of surface tension in 
monovalent metals result directly from the properties of the free electron gas. 
The model used is certainly. fairly adequate for the monovalent metals in the solid 
state, so the surface free energy of the solid at the melting point has been calculated 
for these substances for comparison with the measured surface tension of the 
liquid. The calculated values are about 1-5 times the observed. 


§2. REVIEW OF EARLIER WORK 

The surface energy of metals has been treated theoretically by several workers 
(Frenkel 1917, Gogate and Kothari 1935, Dorfman 1943, Samoilovich 1945, 
Brager and Schuchowitzky 1946a,b). However, only the last two investigations 
are relevant to the present article. | We shall consider briefly the models and 
methods employed by these authors before discussing those of the present work. 

Samoilovich replaced the lattice of ions by a uniform distribution of positive 
charge and employed the '‘homas—Fermi statistical method, according to which 
the total energy of the electrons is completely determined by the electron density 
distribution. ‘lhe positive surface energy thus obtained, however, arises entirely 
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from the use of the Weizsaicker term, which corrects for the kinetic energy 
associated with the rapid change of electron density at the metal surface. This 
may easily be seen, since the original Thomas—Fermi method would give zero 
surface energy for that distribution of electron charge (not a solution of the 
"Thomas—Fermi equation) which is uniform over the region containing the positive 
cloud and zero elsewhere. The distribution which satisfies the Thomas—Fermi 
equation, having lower energy (Lenz 1932), then gives a negative value for the 
Surface energy. 

Apart from the uncertainty of the Weizsicker correction, this model makes no 
allowance for the different ionic fields of different metals, the surface tension being 
determined entirely by the density of positive charge which depends only on the 
valency and density of the metal. 

In view of Samoilovich’s result, that the positive surface energy is due to the 
increase in kinetic energy of the electrons imposed by the steep potential barrier at 
the surface, Brager and Schuchowitzky treated this effect by the direct wave- 
mechanical consideration of free electrons enclosed in a box with impenetrable 
walls, at which the wave function vanishes. ‘They found that the total energy 
contained a term proportional to the surface and identified this term as the surface 
energy. ‘This model neglects the electrical double layer at the surface and, more 
important, the finite values of the work function of real metals. 


Sa, ANSI IMM ONO Nag, 
It is possible to find for the monovalent metals Li, Na, K, Rb, Cs, Cu, Ag, Au 
‘a model which avoids the faults of the two so far considered. ‘Theoretical work 
-on binding energy has shown that in these metals the positions of the conduction 
electrons are so strongly correlated that we may regard each electron as moving at 
any instant only in the field of the ion nearest to it, other ions being neutralized by 
the electrons on them. ‘Then to a good approximation the electrons behave as 
free electrons at a constant potential V, where —eV is the mean potential energy 
between an electron and its nearest ion. ‘The potential naturally vanishes outside 
a metal surface so that we have essentially a Sommerfeld picture of the metal. 
"Thus the model we adopt is as follows. The positive charge of the metal ions is 
‘supposed to be uniformly distributed over a volume proportional to the number of 
ions. Atthe boundary of this volume we imagine a finite discontinuous jump in 
potential. The variation in potential energy of an electron as it passes through 

‘the surface is then as shown in the Figure. ~ 


| 
F | 
Interior of Metal A 


The approximately exponential rise in the potential energy outside the 
boundary arises from the overlap of the wave functions beyond A. There is thus 
a distribution of negative charge in this region which may be calculated by wave 
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mechanics from the electron distribution inside the boundary. The latter is 
calculated by the Thomas—Fermi statistical method. ‘The conditions that an 
equal resultant positive charge must remain inside the metal and that the field 
must vanish in the interior of the metal determine the solution uniquely. 

The electronic surface energy is then divided into two terms : (i) that which 
results from “‘compressing”’ the electron gas within the potential well; this gives. 
a smaller term than one would obtain if ¢ were taken to be infinite as in the Russian 
work; (ii) the electrostatic energy of the double layer : this is found to be small. 


(i) Charge Penetration 


We first find the relation between the charge density outside and inside the 
surface. If we take the x-axis normal to the surface (x=0), the wave function 
representing a stream of electrons with momentum P incident on the surface at an 
angle 6 from the normal and undergoing reflection 1s 


ob, = A(2m¢): exp [1P(y sin 6 + x cos 6+«)/h] 
+ A(2m¢)! exp [¢P(ysin@—xcosA—«)/h], «<0 ...... (1) 
wb, =2AP cos 6 exp [(¢Py sin 6 — Bx)/h], eu 
where tan(«P/h)=f/Pcos@ and B=./(2m¢— P?cos?@). The density p(P, 6) 
of electrons in the incident beam is 
PCP TO) =| Al Zatoe eee (2) 


‘The number of particles per unit area outside the surface is then 
| Iyal2ax=Hip(P, 6)P2cos*8/mgB. eee. (3) 
J 9 


In the Thomas—Fermi method the particle velocities at any point are random 
in direction; thus the density of electrons incident on the surface in the range of 
inclination to the normal from 6 to +d6 and for momentum from P to P+dP is 
given by 

p(P, 0) dP dé =4r7rhP? sin 6 dP dé. 


Thus for the whole Fermi distribution the number 7, of electrons outside the 
surface per unit area is given by substituting for p(P, @) in (3) and integrating over 
P and 6; the result is 


No = 


Z Be [" P*cos26sin@ 
md¢h*) Jo +1/(2md¢— P? cos?) 


-“EAG ams) / Lamta('~ apes) |- G@ -zrg) = (a) 


Here P,, is the maximum momentum in the Thomas—Fermi distribution just 
inside the surface and the dimensionless quantity P?,/2m¢ is the ratio of the 
corresponding maximum kinetic energy of the electrons to the potential energy 


discontinuity ¢. 


(ii) The Thomas—Ferma Equation 


Let E,,, be the maximum electronic energy and V at any point the electrostatic 
potential measured from that in the interior of the metal; then having regard to the 
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constant density of positive charge inside the surface the Thomas—Fermi equation. 
may be set up in the usual way : 


2 
= ae { ean [2m(E,, + eV) 32+ +p Sees 8 (5) 
where 4775/3 =v, is the atomic volume. When the dimensionless quantities 
e=oeeey i,  €=al(fely/3@*)e nea (6) 
where E,, =(97/4)?h?/2mr;, are introduced, (5) becomes 
(Da ae nn ee (7) 


We notice that V-+0 in the interior of the metal, where E,,, = Ey, so that g=1 
| corresponds to the distribution of one electron per atom, as in the interior far from 

the surface. We shall find the deviation of g from 1 to be small in the cases to be 
| treated. It is thus for our purpose more convenient and sufficiently accurate to 


solve (7) approximately to the second order by setting 
eeeiNe Neeser Rie (8) 
and treating A as small compared with unity. 
On substituting (8) in (7) and expanding 93, we find 
Gide A NB Fe ee ee (9) 
when terms of degree higher than the second are neglected. 
The asymptotic solution for large negative € (interior of metal) must clearly be 
of the form 
DATACXD ENS /2)) ASL ee OO. eee (10) 
Thus A is the only arbitrary constant to be determined. The approximate 
solution of (9), to the same order of accuracy as (9), which has the asymptotic 
form (10) is 
A= Aexp [Ev (3/2)] —[A? exp (E/6)]/12 
or g=1—Aexp[&+/(3/2)] + [Az exp(€v/6)]/12. ... (11) 
If the electrostatic field is to vanish at large distances, the sum of the resultant 
charges inside and outside the surface must vanish. This condition determines A. 
The electron density is proportional to g??._ Since when g=1 the density of 
electrons is 3/477}, the net charge density for any value of g is 


3e ae 
eee eee eS ee 12 
o,e Feri g3?) (12) 
and the total charge m,e per unit area inside the surface is thus given by 
3e /? Ap ; 
mem a 0 Eee\dye |) NomaatA, re (13) 
Using (7) and (10), we find on integrating 
me=— pos( Sor) i shouemsieita ( ) 


where gy =(dg/d&)z_0. 

We have now to eliminate P,,, from the expression (4) for the charge outside the 
surface. Just inside, we have P?/2m=E,,+eV=E,gy where go=(g)¢-9, So 
that 7, can be expressed as a function of g : 


mon eB) Lal -8) | - GB) VQ) a 
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where 1/A=E),/¢. It is possible to express gj, gy in terms of A 
&y9=1-A+4+A*/12; oo9= =< |A— Abela), eee (16) 


and A is then determined by the numerical equality of the charges inside and 
outside the surface, 

nA) 0A) ae ee a eee (Tr 
"Thus knowing 7) and ¢ we can calculate the charge distribution near the surface. 
In the cases to be discussed we find that (17) can easily be solved by successive 
approximations, taking as the first, gy=1 in (15). 


$4, LHE ELE CLRICALS DOU Bins “LAYER 

The strength of the surface double layer is measured by the total change in 
electrostatic potential energy of an electron taken through the double layer; it is 
thus given by 

—e[V(0«)—V(—oao)J=—eV(o), — ...... (18) 

‘since we have taken the potential energy to be zero in the interior of the metal. 

To evaluate V(0o), the charge distribution outside the surface may be approxi- 
mated by an exponential function 


ofa) =Cexpie gen ay! peel) tear (19) 


“The density at x=0 is C, which should be equal to the density obtained by inte- 
‘grating the wave mechanical solution | #s.(0)|? as given by (1). Thus we find 


16m (2 \e EPs 
= m m g 
ae (Fi ) Tat sae ae (20) 


‘g is then determined from the total charge : 


n, =| n»| = |. Gish deel, On Wee (21) 


‘The potential V(x) for x >0 can then be determined from the Poisson’s equation 
d?V(x)/dx? =47e0,(x)=4reCexp(—qx).  — .. eee (22) 
At x=(0 the approximate solution must be equal to the potential just within the 
‘surface, given by 
eV) +P 2mak eke or eV,=—f,\ l= 90), eee (23) 


"The solution satisfying this boundary condition and the condition that the field 
vanishes at infinity is 


V =4reC{exp(—qx)-lt/g—-E,y1-—g(O)/e. ...... (24) 
‘So we find for the strength of the double layer 
— CV (co) =4are*C jg? Ha 1 s10)y 5 eee (25) 


We may put this into a more convenient form for calculation by substitution from 


(14), (21) and (22) 
—eV( 0) =Fpgol lt 5go" 20424) eee (26) 


/ . . 
where go, gj are defined as above. Numerical values for the metals considered are 
‘given in Table 1. 
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SO eto es SiURuhy AG Bas BING 4 
The calculation of surface energy on this model must take account of three 
terms: The effect dealt with by Brager and Schuchowitzky (1946a, b) which 
assumes impenetrable boundaries, the penetration of the barrier, and the energy 
of formation of the electrical double layer. The first two terms arise directly from 
wave-mechanical considerations, and are not really physically separable, but it is. 
mathematically convenient to treat them separately. 


(1) The Term of Brager and Schuchowitzky 


This term is due to the increase of kinetic energy consequent on the impenetra- 
bility of the walls. Mathematically this condition is expressed in the requirement 
that the standing electron waves should have nodes at the walls. If we consider a 
cubical box of length L, then in the wave function 4 =sink,xsink,ysink,z the 
wave numbers k,,, k,, k, must have respectively the values /z/L, mz/L, nz/L where 
I, m, n are positive integers. 


Hable 

Metal Cu Ag Au Li Na K Rips 
lB (ieee Taos Veena Ue) D5 5°5 4-7 SOLIS) 2-03 MEATE  al oy 
memeVe yi setiti. 333° 1594 SEZ. 16°3 9-85 8-13 6-48 6:09 5-62. 
S,+S. (erg cm7?) 2180 1400 1490 960 470 208 165 125 
SaiGaeeies) so eel) Mes DD 4:5 iho} 0-32 0-2 EN: 
—F,(ergcm7?).. 360 257 240 66 34 25 Pil 18 
o(ergem~) .. 1820 1140 1250 890 440 180 140 110 
o (exp.) (erg em=2) 1100 800 600-1000 -- 290 200-400 — — 
Double layer (ev.) 0°55 0-29 0-28 0-43 0-25 On2 0:10 0-08: 


iS +S, is the surface energy arising from the change in kinetic energy of the electron gas. 
due to the presence of the surface. 


1s is the surface free energy arising from the change in thermal vibrations. 
S's is the electrostatic energy of the surface double layer. 
o is the surface free energy obtained from these three terms. 


These discrete states are conveniently represented as lattice points in the space 
(k,, k,, k,) with lattice constant 7/L. ‘The number of occupied lattice points, 
multiplied by a factor 2 to allow for the alternative spin states, must be equal to the 
number of electrons to be accommodated and the total kinetic energy is the sum of 


Wee eros 2 IE ae (27) 


over the occupied states. An energy proportional to L? is obtained if we assume 
that the number of lattice points in a spherical octant of volume 7k},/6 is just 
L3k° /7? and that the average energy is simply the integral average of (27) through- 
out the sphere. Clearly this procedure counts also points on the basal planes, 
but with weight 4. These however do not represent permitted states since the 


corresponding wave functions vanish everywhere if they are zero on the walls. 


These rejected points must be made good by taking just half as many at the 
surface of the sphere. The total energy corresponding to those points which were 
wrongly included is 


L\22n (#mh?R? Cie, 
ee 28 
3(= yz Me's =~ kdh= ign ( ) (28) 
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‘The total number of states thus rejected is 3(L/7)*7k?,/4. The energy necessary 
to make them good at the surface of the sphere is 


LP ake Re 
3 (=) mgr hs geet sa trlpel ah  NiaeSe (29) 


Thus the Brager—Schuchowitzky correction, which is the difference of (29) and 
(28), becomes 37(L/7)?h?ki,/16m. Since the total surface of the cube is 6L? this 


m 
is equivalent to a surface energy 


= 7, Clem. pet Util area. a) | |) eee (30) 


Mi 
Brager (1947) has further shown that the same expression is obtained for a spherical 
mass of material, so that its identification as a surface energy is reasonable. 


(11) Term due to Barrier Penetration 


Since the potential energy step at the surface of the metal is actually finite, 
electrons penetrate the barrier and their wave functions are not zero at the boundary. 
In fact each wave function has a definite shift in phase as though it possessed a 
virtual node a little outside the boundary, and there is a corresponding decrease in 
the energy of that state. Let us consider a cube of edge L as in the above section ; 
suppose it to lie entirely in the region x <0 and let the face x =0 alone have a finite 
potential barrier, the others remaining infinite. We see that in equations (1), 
only the factor dependent on xis affected bythe barrier. Ifwewritei-!Pcos@=k, 
this factor becomes 
asin(kyy—) del) “TP = ieee (31) 
where 
tane =k, /./(2imdje— Ry. | eee (32) 


"Thus the wave function hasa phase-shift 6 at the boundary which tends to zero as 
‘b increases. 
The boundary condition at the opposite wall, i.e. at x= —L, requires that 
sin(—k,L —6)=0, whence 
ku= {GO Ee ~~ = eee (33) 


where / is an integer. ‘Thus the substitution of the finite for the infinite barrier 
reduces k,, for each state by 6/L, corresponding to an energy change 


AE =f? A(k2 + k2 + k3)/2m= —(h?/m)k,(8/L). (34) 
"The number of states with R, in the range dR,, is 
doh Rk? Lia) dhe ease Fw eee (35) 
Thus when L is large the energy change per unit surface is given by the integral 
LAE h? Bk k 
S,= => = - 3 hk? —k*)k ,tan-!| —~—* ______ |dx. 
Te nly ee ae a |* 
be mais (36) 
‘The integral is easily evaluated to give 
es, 9 S L 
S,=—- inte {3(2 —A)y/(A— 1) + [3A? — 8A + 8] sin-1 A} 
Pert: (37) 


where, as before, A=¢/Fy. 
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(i111) The Energy of Formation of the Electrical Double Layer 


Apart from the changes in kinetic energy, the energy of electrostatic interaction 
within the system is increased by the formation of the double layer. By some 
authors the surface energy has been attributed entirely to this source. In fact, 
however, we shall see that the electrostatic contribution is entirely negligible 
compared with the kinetic energy terms S, and S,. 

Suppose as before that the surface et in the plane x=0, and let the metal 
occupy the semi-infinite region for which x<0. It is convenient to consider first 
only the finite interval -L<x<L. The energy of the double layer is the sum of 
the self energy of the negative charge cloud, that of the positive cloud, and the 
interaction energy of the two, less the same quantity for a uniform neutral distri- 
bution confined to x <0. 

The distribution of positive charge is unaltered and given by 


PAX) = pes a=): p+(x) =0, hes ie Se eee (38) 


‘Thus its self energy is invariable and need not be considered. ‘The energy of any 
arbitrary charge distribution may be referred to a standard state (zero of energy) 
which may be fixed arbitrarily. We find it most convenient to take the standard 
state to be that in which all the negative charge is concentrated in a surface layer 
atx=0. Then the potential energy between the negative charge in (x,, x, -+dx,) 
and that in (x,—x,+dx,) is 
ATP SP Wo AN ANy| hy Xoj—-— a a tee (39) 

per unit area. For p(x,) dx, is the surface density in the infinite layer of thickness 
dx, at x,, and so 27p(x,) dx,p(x,) dx, is the constant repulsive force exerted on the 
charge in unit area of the layer between x, and x, +dx,, while | x, —x,| gives the 
distance of separation. Similar reasoning shows that the potential energy due to 


the interaction of a layer of positive charge in the interval (x,, x; +dx,) and the 
negative charge in (x, —%*,+dx,) 1s 

— 277% 1p +(%1)p(%2) dx, dx, —2r| xy — X9| p+(%y)p(%z) dx, dx, ...... (40) 
per unit area. Here the first term gives the work done to bring the negative charge 
p(x2) dx, from x =0 to x =x, which is clearly negative. The second term gives the 
work done to separate it again from the positive charge by a distance | x, —x9|, 


which is always positive. The total energy & is obtained by combining (39) and 
(40) and integrating over the various pairs of elements. 


L L 
7 [ax {hs AX {| X, — X9| [ — p(%;)p(%2) — 2p .(%y)p(%2)] — 20 p (%1)e(%2) }- 


For the uniform distribution we have simply 
p(x) = —p,(x). 


Let p’(x) = —p,(x) +Ap(x) be the density of the double layer distribution deter- 
mined above. The energy of the double layer is then 


S3= an oe ie AX) X1 — Xz) 1p(%1)Ap(%2) + 4(%2)Ap(x) + p 4 (%1)p +(%2) 


— 2p .(%,)Ap(x2) — Ap(x1)Ap(%2) = p4.(%1)p (x ‘o) } 
7D cL 
Ee | dee, | daty| x —.v9| Np()Ap(2)). 
25, 
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The second term of (41) drops out directly on forming the difference since the 
total negative charge [p(x)dx, is the same in each of the two cases. Since 
Ap(«) diminishes exponentially in both directions we may pass to the limit L— 20 
and find 
St =| dx, | | diy| ay —¥9| Apa )Ap(m). vee. (42) 

For x>0, Ap(x) is the density of electronic charge which we have called 
—eo,(x). According to (19), it can be represented approximately by the expo- 
nential function —eCexp(—gqx). For «<0, Ap(x) is the net charge density o,e 
defined by (12) which may be written 


Bou, 7" 
~ —-,[3 2 wy), il: tp fone 43 
ae Anr [3 A exp (V/ 2 x)] ( ) 


where we have neglected the termin A?. In the cases considered in the Table the 
inaccuracy thus introduced is of the order of 1% only. ‘Thus we may write 


Ap(x) =O, exp( —qx); - x0; Ap(x)=Oyexp( +G,4), 420 eee (44) 
1676)(Ps eae © Je, 
where OF = = e( h ) 2m¢’ Qy= ee 
= V (3) (Tob p/3e?) rot and Qr/q2=—Qilu «+++ (45) 


The integration in (42) can then be carried out, giving 
61600A? /0:528\3/ 15A ) 
——— ce ene 

A (ion) ( 490" 
where P,, is eliminated using the equation P2/2m=g,F, and V, ry are respectively 
the volume per mole and effective atomic radius in the metal. 

The total electronic contribution to the surface energy of the metal is thus given 
by the sum of the three terms S=S,+S,+S3 where Sj, S:, S; are given respec- 
tively by (30), (37) and (46). 

Only S; involves a quantity which cannot be explicitly given, namely A, which 
depends on the solution of (17), and determines the value of the electrostatic 
_ potential at the boundary, (10), (16). However, we find (see Table 1) that S, 
is only of the order of a few ergs, which is in fact less than the error to be expected 
from a simple model. 


S3= 


(iv) Other Terms 


A contribution to the surface energy from the repulsive interaction of the closed 
electron shells of the metal ions is to be expected. 

This interaction, which is well represented by a central force, would according 
to Shuttleworth contribute about half as much to the surface free energy per atom 
as to the binding energy. ‘hus if it contributes a rather small fraction to the 
cohesion, it will in the monovalent metals contribute between three and four times 
that fraction of the surface energy. This is negligible both for the bulky alkali - 
metals and for the noble metals with their hard completed d-shells (Fuchs 1936, 
Huntington and Seitz 1942). 

A more important effect is the change in character of the lattice vibrations of the 
metal in the presence of a free surface.  Brager and Schuchowitzky (1946 b) have 
considered an extension of the Debye model of a crystal to a large thin sheet of 
material with free boundaries, and obtain a term proportional to the total surface 
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area in the thermal energy. This surface internal energy due to the lattice 
vibrations is given by them as 


2: 15 kT 4 u-du x3 2 
Ve Xe [2| pees =| ergs/cm*, == wee eee (47) 


= 


where x=©/T and © is the Debye characteristic temperature of the material. 
The corresponding value of the free energy is 


°T 
F,=E,-T | T-\(dE./dT) dT 
~ 0 


which can be reduced to 
Ne Soe i VACC} 4) al ene lees Pre (48) 


where A,, =2-15RO/v?? and Z(x) =x an wee 48° that A,, is a constant for each 


_ metal determined by the Debye characteristic temperature and the atomic volume. 


The integral occurring in the function Z has previously been tabulated 
(Zanstra 1931), but only for rather large intervals of the argument. Since Z(x) 
varies as 1/x for small x, it was decided to recalculate the integral and Z. The 


results are shown in Table 2. 


x 2 
@ablew2 sila) — = = ; Z(x) =x-8 I(x) 
/ 0 bag 
ee I(x) Z(x) oy I(x) Z(x) x I(x) Z(x) 
0-1 0-004 835 4-835 0-4 0:069 865 1-092 0-7 0-192 808 0-562 
0-12 0-006 916 4-003 0-42 0-076 498 1-033 O72 5 0-202)5585 10-543 
0-14 0-009 351 3-408 0-44 0-083 381 0-979 0-74 0-212 471 . 0-524 
aon 0-012, 131 2-962 0-46 0:090 507 0-930 O27 Or 02227 543es URS OV) 
Orton 2 0-015,250 > 2°615 0-48 0-097 871 0-885 O77 Se O2232) 167. OFF 
| 0:2 0-018 700 2-338 0:5 0-105 465 0-844 0:8 0-243 139 90-475 
0:22 0-022 474 2-111 0-52 0-113 283 0-806 0:82 0-253 654 0-460 
0:24 0-026565 1-922 O:54 0-121 321, 20-771 0:84 0:264 307 0-446 
0-2) 0-030 966 1-762 0-562 201295 72a— 07358 0:86 0:275 093 0-433 
0:28 0:035670 1-625 0°58 0-138 030 0-707 0:88 0-286 008 0-420 
0:3. 0-040 668 1:506 0-6 .0:146 689 0-679 0:9 0:297046 0-408 
0-32 0-045 956 1-403 0:62 © 0°155 543. -0:653 0-92 0:308 204 0-396 
0:34 0-051 527 1-311 0-64 0:164 588 0-628 0:94. 0-319 477 0°385 
0:36 0-057 373 1-230 0-665 0-173) 81-7, 20605 0:96 0:330859 0-374 
0:38 0-063 488 1-157 0:68 0:183 226 0-583 0:98 0-342 348 0-364 


It should be noted that the contribution to the total surface free energy from 
this source is negative. 

This calculation of the surface free energy due to lattice vibrations is, of course, 
only valid for the crystalline solid. For the liquid, consideration of capillary waves. 
(Frenkel 1946) provides the most reasonable approach, but in the absence of a 
more detailed theory of liquid structure such considerations cannot be subjected 
to quantitative test. 

We have therefore calculated, for approximate comparison with the experi- 
mental results for the surface tension of the liquid metals near their melting points, 
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the total surface free energies of the solid metals at their melting points, using the 
expression 


SENS ES 


“IE Lfsa oven sen-sense (CY) 


=2-11(0/V**)Z(0/T,,) ergsjem™ 0) 0 are -. (49) 


where V is the mole volume in cubic centimetres and ©, 7, the Debye characteristic 
temperature and melting point respectively in degrees absolute. A is the ratio of 
the barrier height ¢ to the maximum Fermi energy Ep. 


§6. NUMERICAL RESULTS 


Theoretical work on the noble and alkali metals has shown that the conduction 
electrons of these metals behave almost like free electrons. ‘Thus we can use the 
ordinary electron mass in our calculation, so that the electronic contribution to the 
surface energy depends only on the height of the potential barrier ¢ and the mole 
volume V. The estimate of ¢ most consonant with the model we employ is obtained 
from the consideration that the binding energy W of the metal (atomic ionization 
energy plus heat of sublimation) is related to ¢ by the equation (Huang 1948) 


-W=3(-64 7) SNC ie oi ae (50) 


The calculated values of the various terms in the surface energy and the strength of 
the double layer for the noble and alkali metals are given in Table 1, together with 
the experimental values used in the calculation. ‘The experimental values of the 
surface energies of the liquid metals as given in Brager and Schuchowitzky (1946) 
are given in the Table for comparison. 

The agreement with the theoretical values is within a factor of about 1-5, 
which is as much as can reasonably be expected from this model. The relative 
values for different metals agree well. ‘There are no definite experimental values 
for the double layer moments, but the values found are in substantial agreement 
with more elaborate theoretical calculations obtained (Bardeen 1936) in this theory 
of the work function. 

One might notice that the differences (f — E,) formed from the values of ¢ and 
E, given in Table 1 are much higher than the experimental values for the work 
functions of the same metals. ‘This discrepancy is, however, only apparent. It 
is a consequence of the fact that we have used a single barrier height ¢ for all 
electrons, whereas in reality different electrons have effectively different barriers 
owing to the difference in the exchange correlation energy —«(k) for electrons with 
different wave number k (cf. Seitz 1940). Thus the effective barrier for an electron 
with wave number k is higher than that for an electron at the top of the Fermi 
distribution (wave number kp), by <(k) —e(ky). If x is the observed work function, 
so that the barrier for an electron at the top of the distribution is ¢(ky) = x + Ey, the 
effective barrier for an electron with wave number k is 4(k) =y + Ey is — (Rp). 
The mean value for all electrons is 


$=xXt+Ey +e —e(ko) =x + Ly +0-916€7/7)—-0-712e2/7, ...... (51) 


where the expressions for €, «(ky) are from Seitz (1940). 
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If the mean barrier heights thus determined are used in the calculation, results 
for the surface energy are obtained which differ only very slightly from those given 
in Table 1. 
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ABSTRACT. Observations on the spectrum of sodium hydride have been extended to 
longer wavelengths (from 4600 a. to 6450 a.) and photographed under high dispersion. 
Analysis of the data thus obtained gives values of the vibrational and rotational constants for 
several new vibrational levels, among which the lowest levels of the excited-state are of 


especial interest. 

The intensity distribution has been calculated from the wave functions appropriate to 
potential energy curves based on the constants derived from the analysis of the band system. 
"Tne theoretical estimates are found to be in fair agreement with the experimental results. 


§1. INTRODUCTION 

HE absorption spectrum of sodium hydride was reported by Hori in 

| 1930. The values assigned to v’ were subsequently raised by one unit 
as a result of work on the emission spectrum (Hori 1931), and later by a 
further three units on the basis of a comparison with the deuteride (Olsson 1934). 
No data were available for these lowest three vibrational levels of the excited 
state. They are of especial interest because the values of w, and BY’ for the 
hydrides of the alkali metals behave anomalously in that they first increase with 


* Now at the National Physical Laboratory, Teddington. 
14-2 


192 R. C. Pankhurst 


v' instead of decreasing throughout. The present work was undertaken with a 
view to determining the missing data; the preliminary results have been reported 
elsewhere (Pankhurst 1941). 


§2. EXPERIMENTAL 


The spectrum was obtained from a hydrogen discharge tube source of the 
type described by Pearse and Gaydon (1938), the discharge passing through a 
quartz constriction containing a pellet of sodium. The tube was cooled by 
immersion in water, a current of from 0-3 to 1:0 amp. being supplied from a 
2,500 v. transformer. , 

The band system was photographed with 6-m. gratings, using Ilford H.P.2 
and hypersensitive panchromatic plates. Between 4600 a. and 5660 a. a 
dispersion of 1-9 a/mm. was obtained by means of the second order of a grating 
with 11,000 lines/inch (433 per mm.); further to the red the first order of a 
grating with 15,000 lines/inch (591 per mm.) was employed, giving a dispersion of 
2-6 a/mm. The whole system was also photographed on a Hilger E2 quartz 
prism spectrograph. 


§3. ANALYSIS OF EXPERIMENTAL RESULTS 


The system is of the many-line type, with no apparent band-heads. It 
arises from a !X—'> transition and resolves, on analysis, into R— and P— branches. 
Some 25 new bands were identified, involving transitions from seven vibrational 
levels. (A copy of the detailed rotational analysis of each band may be seen 
on application to the author.) 

Rotational and Vibrational Analysis. Mean values of the observed combination 
differences, obtained by averaging the values for all the bands obtained from each 
of the vibrational energy-levels, are set out in Table 1. The rotational analysis 
was based on the usual power-series representation of the rotational energy, 
the constants thus derived being given at the foot of the Table. The closeness 
to which the experimental results are represented is shown by the agreement 
between the observed combination differences and those calculated from the 
derived rotational constants: the agreement is everywhere well within the probable 
error of the observations. 

The vibrational analysis is set out in Table 2, which includes Hori’s (1931) 
values for bands involving transitions to the lowest three levels of the ground 
state. The vibrational energy in each of the two electronic states is represented 
by the power-series expressions for G(v+ 4) given in the Table. The accuracy 
to which these expressions fit the experimental data is indicated by the agreement 
of the calculated band-origins included in Table 2. 

Note on Accuracy. It is important not to claim high accuracy for the absolute 
values of constants obtained as coefficients in power-series expansions such as 
those for F(K), By(u), D,(u), H,(u) and w,(u), because a substantial alteration 
to any one of the coefhcients can be largely compensated for by appropriate 
adjustments of the others. It is in fact possible to find several polynomials of 
the same degree but with corresponding coefficients considerably different, 
both of which represent the original observations to within the experimental 
error. Although this may be immaterial if a power series is to be used merely - 
for the purpose of interpolation, it introduces considerable uncertainty when 
(as in this present case) the coefficients are to be interpreted as physical properties. 
of the system. . 
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Table 1. Rotational Analysis 
(a) Excited State 
(1) Mean values of the observed combination differences, Aj(K) = R(K)— P(K) 


y 
ae 1 2 3 4 5 6 | 
K 
1 11-85 10-86 11-52 11-69 12-23 
2 18-46 18-68 19-42 19-33 19-77 19-38 
3 25:22 26:16 26:87 26-94 27:43 27-47 
4 32-66 33-61 34:20 34-57 34-84 34:89 34-68 
5 39-47 40-91 41-60 42-17 42-43 42:25 42-37 
6 46-84 48-26 49-13 49-68 49-95 50-06 50-02 
7 54-00 55-54 56-46 57-15 57-49 57:55 57:43 
8 60-94 62-67 63:82 64:57 64-79 64-99 64-91 
9 67-76 69°72 71-06 71:81 72:37 72:33 72-25 
10 74-54 76:70 78-17 79:10 79°55 79-68 79-56 
11 81-21 83-59 85-19 86-14 86-60 86-87 86-72 
12 87-82 90-31 92-07 93-16 93-66 93-94 93-87 
13 94-23 96-94 98:76 100-03 100-70 100-96 — 100-85 
14 100-63 - 103-46 105-37. 10680 107-49 107-74 107-67 
15 106-88 109-83 11200 113-41 114-17 114-49 114-38 
16 113-00 11618 118-38 119-90 120-68 120-84 121-06 
17 fis 07 eto? 19 e467" 2126-29 12713. 127-51. | 127-56 
18 124-84 128-22. 130-67 132-46 = 133-50: 133-95. 133-83 
19 130-65 13410 13669 138-44  139:50 140-03 139-97 
20 136-42 139-76 142-47 144-24 145-42 145-98 145-97 
21 141-75 145-40 148-11 150-210 451229 9 4'51-96" 154-78 
22 ay 1oe 150-75 1 453-66" 155-6372 -156-97  - 158-18 157-42 
23 t2-O8_ 9156-16. 15912 1161-15 = 7 162-38 162-79 
24 157-54 161-37. 164-33. «166-42: 167-71 168-26 
25 162-70 166°37 169-31 171-49 172-73 172-95 
26 167-60 171-29 174-24 17646 178-01 
27 176-13. 179-06 181-36 182-70 
28 180°80 = 183-73 187-50 
29 185-32 188-20 191-92 
30 189-33 192-61 196-10 
31 196-90 


(ii) Values derived for the 


v’ 1 2) 3 

1B 1-823 1:875 1:908 
Dix 104 .| —2°25 —2:25 —2:20 
ee Se 1108 5-10 3°85 3-04 


rotational constants (in cm!) 


4 5 6 7 

1-930 1-938 1-941 1-936 
= 29 oa (i7 2-0 193} 

2°87 2-02 Led 1-24 


Fi(K)=B(K(K+1)+D\K%(K+1)?4+ Wy KK -+ 1). 


The above constants may be represented by the equations : 


Re Bad pag? 6... (u’=0’ +4) 


‘ 2 
D,=Do-+ Bu’ + Baw? ...- 

[ey oo ff a, 7,,72 
He=HHe- yu tyr fe. 


with the following values of the coefficients 


Bi=1-696, whence J=16:31X10-gm.cm? and 7¢=3-:200a.; a; =0-1083 ; 
a,=—0-0175 ; a,=0-00129 ; a,;=0°000042 ; Dg=—2:27x10-* ; i= —0:01 x 10-4 ; 
e008 104 = H5=5:7X102; y,;=—0°65x10-. 
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Table 1 (cont.) 
(b) Ground State 


(i) Mean values of the observed combination differences, 
A3(K) = R(K — 1) — P(K +1) 


: aie 
3 a 5 6 7 8 
vom : 
1 26°45 24-90 Dey | 
2 44°33 42°50 42:07 41-33 39-27), 
3 61-98 60°22 58-39 56-70 54-70 
4 LOY TT 24 74-96 72°58 70°46 
5 9713 94-32 91-48 88-60 85-83 83-01 
6 114-63 ters 107-98 104-62 101-14 97-94 
7 13196 128-06 124230 120-34 116-59 112-66 
8 9-74 144-76 140-45 136-12 131-74 127-30 
9 166-41 161-37 156-53 151-83 146-94 144-77; 
10 183-27 7-84 172-50 167-16 161-77 156:45 
11 200-10 HOA 188-27 182-41 176-62 170-69 
i, 216-73 210-35 203-84 197-49 19127 184-73 
13 233.214 226:23 21932 212-46 205-70 198-74 
14 249-40 241-93 234-63 227-20 219-84 212-52 
5 265-49 D5 5y/ 249-68 241-78 233293 226-05 
16 281-31 272°81 264-49 256-10 247-77 239-36 
17 296-76 288-02 DU sV7 270-19 261-34 252-43, 
18 Size) 302°78 293-36 284-03 274-71 265-35 
19 327-24 317-26 307-60 29757, 287°77 211-99: 
20 342-03 331-61 321-24 310-93 300-62 290-22 
72) 356-41 345-52 Saye 7) — SUS SU57 
22 370-70 359-40 348-10 336-84 3257335 
23 384-61 SMe? 360-98 349°45 335i 
24 398-25 385-92 SA SaEA 361-15 
2» 411-55 398-68 386-02 
26 424-45 411-13 398-07 
DY 437-15 423-30 409-65 
28 449-32 AS aeil() 421-00 
29 461-10 447-28 432-04 
30 472:63 = 442:-63 
31 483-79 468-72 452-75 
(11) Values derived for the rotational constants (in cm~!) 
vo" 3 4 5 6 7 8 
iBe 4-438 4-302 4-180 4-051 SNS 3-790 
Dy x 104 —s.3il Os = 3°29 Sel = S05 — 3°04 
H,, x 108 Zl, desl 2:58 2°70 0-69 1:10 


Fi(K)=B'K(K+1)+D'KXK+1)?-+ H¢K(K-+1)3. 
The variation of By may be represented by the equation BJ =Bg+aju", (u’=v" +4), 
with Be=4-886, whence I¢=5:6610-4° gm.cm? and 7(=1:885a.; aj=—0-129. 
Mean value of Dy=—3-15 <x 10-4 ; mean value of H?=1-7x 10-8. 
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Energy of Dissociation. The usual equation 
DA eA ee (1) 


assumes a regular closing up of the vibrational energy levels. For NaH this 
condition is satisfied in the ground state but not in the excited state, for which 
«, at first increases with v. The above equation was therefore used only for the 
ground state; the value for the excited state was obtained from that for the ground 
state by assuming that the products of dissociation were a normal sodium atom 
and a normal hydrogen atom from the ground state, and an excited (?P) sodium 
atom and a normal hydrogen atom from the excited state. The results were: 


D’ =11,650 cm-!=1-44ev.; D”’=17,410 cm! =2:16 ev. 


§4. POTENTIAL ENERGY CURVES, WAVE FUNCTIONS AND 
INTENSITY DISTRIBUTION 


In order to estimate the intensity distribution by the methods of wave mechanics 
it is first necessary to determine the potential energy curves for the two electronic 
states and then to calculate the wave functions for all the vibrational levels. 
If % and %” denote wave functions for the excited state and the ground state 
respectively, the intensity of each band (v’, v”) is obtained from the expression 
f(W’b" + h'b") dr, where the bars denote complex conjugates. 

Excited State. ‘The potential energy in the excited state was written in the form 


VAG ae earns ane Cee. ae le Se SR (2) 
where 

Fi=(7—7,) in, Se eee oe (3) 
with the coefficients evaluated from the vibrational and rotational constants 
using the expressions given by Dunham (1932) but neglecting the small correction 
terms of the order B?/.?. The values obtained were: 


mae a! = 0-949 a= 3:28 ah = 7-09. 


The resulting potential curve lay entirely within the parabola ajé’? and 
resembled another parabola of smaller latus rectum. This observation was 
exploited when calculating the wave functions ’ (§5). For most molecules 
the right-hand limb of the potential curve falls below the simple parabola; the 
irregularity in this instance is associated with the anomalous behaviour of the 
constants for the excited state. 

Ground State. Inthe ground state the potential curve given by a similar power 
series flattens prematurely, for positive values of €, between v” =4 and v’=5. 
As dissociation could not have taken place at these values of v", it was concluded 
that an expression of the form cf equation (2) is not suitable for representing the 
potential energy of NaH in the ground state. 

The expression due to Morse (1929) was therefore used instead: 


: | Vani ltane remy Se Bones (4). 
with D” = 17,410 cm™ and a” =1-:070 at. 
Intensity Distribution. ‘The intensity of the band v’>v" is 
ae 5 
I(v', v")=C exp {— heG'(v')/RT} 4 | (cb b" + p's") ar | Bee Se (5) 

0 

where C is an arbitrary constant, k the gas constant per molecule and 7 the absolute 

temperature. 
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§5. COMPUTATIONAL DETAILS 
Excited State. Asan analytical solution was not available for the wave functions 
appropriate to a power-series potential curve, they were obtained for the excited 


state by modifying those for the parabola ajé’, for which the solution may be 
written 7 =bF(v, r) where 6? =27/(wew,/h) and 


Po )= 1 (sae )exw(-8 2) 


uw being the reduced mass of the molecule, and H, denoting successive Hermitian 
polynomials with s =b?(r —r,). 

The modification consisted simply of adjusting the abscissae so as to “ fit’” 
the wave functions to the limbs of the power-series curve instead of to the limbs 
of the parabola. ‘The ordinates for each vibrational level were then multiplied 
by aconstant factor in order to maintain normalization.’ ‘This procedure, although 
somewhat crude, is to some extent justified by the resemblance of the power- 
series potential curve to another parabola of smaller latus rectum (as already noted 
in §4), and might be expected to give results at least as accurate as the eye-estimates 
of the observed intensity distribution. As a check, a parabola was fitted to the 
power-series curve by the method of least squares ; its wave function was evaluated 
for one vibrational level (v' =5) and was found to agree reasonably well with that 
obtained by the arbitrary method indicated above. 

Ground State. The wave functions corresponding to equation (4) have been 
derived by Morse (1929) and may be written in the form 


Py ee Pes ee ig. ee eee (6) 
where Po(hay et SS ae (7) 
T(k=v 
Rate exp {im(k —v—1)}, Ainicac. S45 (8) 
Fp =e 2a(k-2 D2 ye — (ky — 12" + (k--v—1)R > - —2)o(v—1) ge2— | 
roeuee (9) 
and k=4n(2uD)?/ch; S=—RCXPY=-aAr=7.\y. | Beaeeee 10) 


N,y=[P(R-v)]? & D(R-20+p—1)/T(p-1). 
p=0 

The numerical computation of equations (6) to (11) would involve a prohibitive 
amount of labour, partly owing to the large number of points needed to define 
the wave curves for the higher vibrational levels, and partly because each term in 
equation (9) has to be evaluated to extremely high accuracy in order to obtain 
relatively few significant figures in the sum, since the positive and negative terms 
are so nearly equal and opposite. Equation (6) was therefore evaluated for only 
three of the levels (namely v” =0, 4 and 6) and an approximate method was devised 
for circumventing the numerical evaluation of equation (6) for the other vibrational 
levels. 

This approximate method is suggested by the observation that if the expo- 
nential in equation (4) is expanded in series, the first term is of the form aj€2, 
for which the wave functions can more readily be obtained. For the remaining 
levels of the ground state, therefore, wave functions were first calculated for the. 
stmple parabola and then modified in a manner similar to that adopted for the 


Pee oe ee 
i ee ee 
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excited state. Normalization was preserved by scaling the ordinates by the con- 
stant factor appropriate to each vibrational level. The method was justified 
by a comparison of the values of (f’" +s") (cf. equation (12) below) with those: 
obtained from the exact solutions for v” =0, 4 and 6. 

_ Intensity Distribution. With the potential energy curves of equations (2) and (4), 


| Bp" +p yp" =b'y", STi ae R (12) 
where’ y"= FU F, Fj y.; with Ff and Fé , given by (7) and (9) but with. 
" ” eal r k” =o" " wu 
Few LN yy — vacate it — Lanai ake Aecstaars (8a) 


In evaluating equation (5), a temperature of 680°c. was assumed, corre- 
sponding to a value of 1-5 for he/RT. This assumed temperature may be a. 
little high, but will not materially affect the general comparison with the experi- 
mental results. The arbitrary constant C was so chosen as to give values of 
intensity on a scale of 100 for the strongest band. 


$6. THEORETICAL RESULTS 

The calculated intensity distribution is compared with the experimental. 
eye-estimates in Table 3 (inset in Table 2). The theoretical values show the three 
observed loci of maximum intensity indicated by heavy type in the Deslandres. 
array. Moreover, the positions of these loci are given with fair accuracy, despite 
the high values of v’ and v” at which the secondary maxima occur, where the 
approximate wave functions used might be expected to be least reliable. In 
view of the several approximations introduced in the numerical work, and the 
various assumptions involved in the basic theory, the general agreement with the: 
observed intensities is surprisingly good. 

The trend of the primary (Franck—Condon) locus of maximum intensity 
for small values of v’ explains completely the drastic extension of the spectrum 
which had proved to be necessary in order to obtain bands involving transitions. 
from the lower vibrational levels of the excited state. 
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Bending Vibrations of a Linear Chain 
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ABSTRACT. The bending frequencies of a linear chain of m mass points have been 
calculated using a harmonic potential function involving forces between second neighbours. 
The proper boundary conditions are introduced and a general formula for the bending 
frequencies is derived. 


H longitudinal vibrations in a linear chain consisting of equidistant mass 

points held together by quasi-elastic forces have been known for some 

time although no attempts have been made to compute the bending fre- 
quencies of a linear chain containing a finite number of mass points. In the 
case of a continuous distribution of mass points, this case passes into the case 
of the transverse oscillations of an infinitely thin rod. 

The linear chain will be assumed to consist of m equidistant mass points 
of mass m which are held together by valence and angle forces, i.e. forces between 
immediate and second neighbours. The potential function V, assumed to be 
harmonic, of such a chain will then be 

ew ik 


Vier 
DV afi 1s ade SA) bie 9 SY eee ee (1) 
==] i=2 D 


‘where the J; are the variations in the distances between P,,, and P,, the 4; are 
the angles between P,_,, P; and P,,, produced by the oscillations, and a is the 
equilibrium distance between the mass points. If u;, v; are the horizontal and 
vertical displacements respectively of the vibrating mass points, it 1s easily seen 
that 
[= Uy — 3 A= (—1)K20-Ga 4)! ee (2) 

assuming that u;2 and v,? are very much less than a?. It follows that the vibrations 
of such a chain split up into two branches, that of the valence vibrations due to 
the force f between adjacent masses, and that of the bending frequencies produced 
by the force ¢ between second neighbours. The valence frequencies arise from 
the longitudinal displacements of the particles and the bending frequencies 
from their movements normal to the direction of the chain, in strict analogy 
to the longitudinal and transverse oscillations of a bar or rod. 

The differential equations of the transverse vibrations are, according to (1) 
and (2) 

— m0; = (—1)'ab(2A,+A_3 +241) 
or — Md; = $(Y_2 —40;_1 + 69; — 40,43 + U4): 
If we write A@v,=v;,,-u,; A®v,=9,,.. —2u,,, +9; and in general 
k a, — > ya) 
AX VO NN DO WE Dy, 

the system of equations (3) can be written in the form of difference equations: 
and ai becomes (= AS ae Equation (4) is analogous to the differential 
equation determining the transverse vibrations of a bar py = —xq(d4y/dx") where 
p is the density of the material, y the vertical displacement of the point with the 
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abscissa x, « the radius of gyration about a line through the axis, and q Young’s 
modulus (Rayleigh 1929). Using the expression (2), the system (4) transforms into 
st pl Nga 4A, OAs a Agi Ata), ea ede: (5) 


There are n —2 simultaneous equations corresponding to the m — 2 valence angles. 
The boundary conditions can now be derived from the equations for the two- 
initial and final angles which have to be evaluated separately from (1) and (2). 
These equations will be somewhat different from (5) since the points P,, Po, 
P,,_, and P,, have not the full complement of second neighbours. ‘There result. 
the equations (6): 
— ms = (6s + 43 +A4) 5 — mh; = $(4A, + 6A3+ 4A, +A5); 


f MA, 9 = P(An—4 + 4X3 a 6A,» a aes} ) aa mr, 4 i P(An-3 a6 4n, 9 oa GN, sa) i 


These four equations can be said to have the same form as (5) if it is assumed that 
Xo; Ay, Xny Ap 41 =9 ) 
or NOG, AOU, A, 7, AO, = 0.) 
The equations (6a) represent the four boundary conditions ; they indicate that no- 
force couples act on the two ends of the chain. In the case of thecontinuous 
bar they reduce to the conditions that d*y/dx*=0 for x=0 and x=1 if the two 
ends of the bar represented by x =0 and x=/are free. A, is now assumed to be 
d, = cos wt{Ae™ + Be” + (—1)'Ce"? +(—1)"De*}._ If the first three boundary 
conditions are applied, the following value obtains for A;: 
\;=A’' cos wt{(—1)"*' sine sinh (n —2)8 —sin (nm —2)x sinh B 
+(—1)" sinia sinh (m—1)6 +(—1)”" sin(—1)« sinh »B 
—(—1)'sinn« sinh (z—1)8 —(—1)'sin(m—1)x sinhzB}. ...... (7) 
The last boundary condition leads to 
0=2 sinw sinhB +(—1)" sin(m+1)a sinh (m—1)6 
+(—1)” sin(m—1)« sinh(m+1)8+2(—1)" sinna sinhnB. ....... (8) 
The relation between « and f can be determined in the following way: if equation 
(5) is solved by a term Aa’ cos wt for A,, it follows that a is a solution of the equation 


a+4a*+(6¢6—mw*)a/ot+4at+l=0. —...... (9) 
If the roots of (9) are aj, dy, a, and a,, we have 
4 4 4 
4 a= —4, G1, La;= LD 4,a,a,,. 
i=1 i=1 i=1 iAj#k 
hese relations are fulfilled if a,=1/a,, a,=1/a,. If then a,=e%, a,=e-%, it 
follows that a;+1/a,= —4—2cos%. If a; is assumed to be —e’, the relation 
between « and f is 
COSMO =a COSneT ™ ae (10) 
Expression (9) shows that 66—mw?/$= &% aa,=2(1—2cos« coshf), therefore 
i#j 


w =27yv can be determined, taking (10) into account: 
mw? =4¢(1+cosa)?; v=[(1+cosa)/m]./(d/m) O<a<m. ....... (11) 


For every value of n, there are 2 —2 solutions for cos«. All the possible bending: 
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frequencies are contained within the limits v,,,,=0 and »,,,.=(2/7)\V/(¢/m). 
£ can-assume values between 0 and 1:76. The bending frequencies of a linear 
chain are doubly degenerate in the two planes through the axis so that the 7 —2 
frequencies (11) represent 2n—4 vibrations. The n—1 valence vibrations are 
given by the well known formula 


int, 
yah ff, oak ny REV 28a ode Ne Gee (12) 
7 Vit 


The 3n—5 possible vibrations of a linear chain consisting of m mass points are 
therefore represented by the equations (11) and (12) provided that the particles 
are subject to valence and angle forces. 
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Sessile Dislocations 


Let us consider what would at first sight seem to be the simplest sort of dislocation 
in a close-packed cubic lattice, made, in imagination, by cleaving the crystal part way between 
two close-packed (i.e. {111}) planes and inserting a close-packed half-plane into the cleft. 
‘This can be done with maintenance of close-packing of nearest neighbours for all atoms 
except those at the edge of the half-plane, i.e. along the dislocation line itself. This would 
appear to be an ordinary Taylor dislocation (Taylor 1934). However, should this disloca- 
tion now glide, in the ordinary manner of a Taylor dislocation, in the direction normal 
to the added plane, it will bring an ever-increasing number of atoms out of the close- 
packing arrangement into another arrangement of more open packing and consequently 
of higher energy: so that this glide is prevented by a large restoring force—so large, we 
may easily conclude, that the glide would not take place under any macroscopically applicable 
stress. Such a dislocation will be called “‘ sessile’, in contrast with “‘ glissile ’ dislocations— 
those which are capable of glide. 

The situation may be clarified by considering the displacement vectors—Burgers 
vectors (Burgers 1939)—of various possible dislocations. The displacement vector of a 
perfect dislocation is a lattice vector, ($a, $a, 0), with its various permutations, in the 
close-packed cubic lattice. Besides these perfect dislocations, which transfer atoms 
from one lattice position to another, there are also imperfect dislocations which transfer 
atoms from lattice positions. to twin-lattice positions. Such are Shockley’s “ half- 
dislocations ’’ (Heidenreich and Shockley 1948) with Burgers vector (a/3, a/6, a/6) which 
arise by the dissociation of perfect dislocations in a {111} plane, e.g 


sa, —yza, 0)—>(a/3, —a/6, —a/6)+(a/6, —a/3, a/6). 
Such is also the dislocation with Burgers vector (a/3, a/3, a/3) described in the first para- 
graph of this letter, and one with Burgers vector (—a/3, —a/3, —a/3) obtained by removing 
instead of adding a close-packed half-plane. Unlike perfect dislocations of opposite 
sign these two are not equivalent to each other. 

Now, atoms in twin-lattice positions only fit correctly on the appropriate composition 
plane (a {111} plane in these cases) of the main lattice. Hence, whereas a perfect dislocation 
can glide, without changing the amount of misfit in the crystal, in whatever surface contains 
the line of the dislocation and its Burgers vector, imperfect dislocations increase the misfit 
unless they glide in a composition plane. Consequently the only kind of imperfect disloca- 
tion which can glide freely is one whose Burgers vector lies in the composition plane of 
the twin relationship to which it belongs. The imperfect dislocations of (a/3, a/6, a/6) 
type, lying in {111} planes, are glissile. The imperfect dislocations (a/3, a/3, a/3), not 
lying in {111} planes, and also being incapable of dissociating into glissile partial dislocations, 
are sessile. The only simple way in which they can move in the lattice is by diffusion 
processes in which they gain or lose atomic vacancies. 
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The ways in which these sessile dislocations can be formed (e.g. by the aggregation 
of vacancies, or by relatively complicated dislocation reactions) will be the subject of a 
longer communication. It is thought, however, worth while to draw immediate attention 
to their existence. ‘They have obvious importance in connection with the strength of metals 
(being fixed sources of stress which exert forces on glissile dislocations) and with crystal 
growth (Burton, Cabrera and Frank 1949). The sessile Taylor dislocations here described 
produce the same type of permanent terrace on the {111} faces they intersect as is produced 
by an (a, 0, 0) screw dislocation in the {100} face of a simple cubic crystal model. 


H. H. Wills Physical Laboratory, F. C. FRANK. 
University of Bristol. 
22nd December 1948. 
Burcers, J. M., 1939, Proc. Kon. Ned. Akad. Wet., 42, 293. 
Burton, W. K., Casrera, N., and FRANK, F. C., 1949, Nature, Lond. (in the press). 
HEIDENREICH, R. D., and SHOCKLEY, W., 1948, The Strength of Solids (London: Physica! Society), 
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Taytor, G. J., 1934, Proc. Roy. Soc. A, 145, 362. 
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ABSTRACTS FOR SECTION B 


A Photoelectric Type of Acoustic Spectrograph using Sound Film, by D. Brown, 
C. F. CoLteman and J. W. LYTTLETON. 


ABSTRACT. An account is given of an apparatus which analyses non-recurrent wave- 
forms such as occur in speech, so that frequency-distribution is represented against time 
in an FT-diagram. From the theory governing such diagrams, no unique representation 
is possible, the result depending on the analysing convention adopted, but all conventions 
are governed by the uncertainty relation AFATS1. The method employed is an optical 
one which analyses a continuously moving sound film by means of a frequency-scanning 
disc, so producing an oscillograph pattern which is photographed on to another moving | 
film. Typical analyses are reproduced and discussed, mainly “ speech patterns”? which 
show the composition very clearly. Advantages include the ability to record the 
amplitude and phase of each component, also the ability to change the analysing convention 
in any desired manner, for example to give the error-function convention which leads to 
minimum uncertainty. The paper concludes with considerations of resolving power, 
spectral line profiles, and amplifier band width requirements. ; 
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The Cause of Anisotropy in Permanent Magnet Alloys, by K. Hosetitz and — 
M. McCaie. 


ABSTRACT. Measurements of magnetostriction have been made on samples of an 
anisotropic permanent magnet alloy. From the results it is concluded that the domain 
magnetization is, in the absence of a field, along that easy crystallographic direction which 
makes the smallest angle with the axis of anisotropy. This view is supported by measure- 
ments of the remanence. The ideal arrangement which obtains after hardening seems 
to be slightly disturbed during the subsequent tempering operation. 


Magnetic Behaviour of Zinc and Cadmium at Low Temperatures, by L. MacKINNON. 


ABSTRACT. Following the recent discovery of the de Haas—van Alphen effect in zinc 
by Marcus at liquid hydrogen temperatures, the effect has been investigated at liquid helium 
temperatures using (for the most part) a torsion method similar to that used by Shoenberg 
in investigating bismuth. The effect in zinc takes the form of a variation of the susceptibility 
along the hexagonal axis which is periodic in the applied magnetic field. An attempt has 
been made to correlate the results with Landau’s explicit formulae and to derive the various 
relevant parameters (such as effective masses, and degeneracy temperature) for the electrons 
causing the effect. Only qualitative agreement has been found, so that too much reliance 
should not be placed on the physical significance of the figures derived; however, it would 
appear that, as in the case of bismuth, the effect is caused by only a few times 10~* electrons 
per atom. In spite of its very similar structure, cadmium did not show the effect; it had, 


however, at liquid helium temperatures, a magnetic anisotropy some six times that at room 
temperature. 


On the Dielectric Properties of a Gas Discharge, by E. E. eS and R. E. B. 
MAKINSON. 


ABSTRACT. Formal expressions are derived for the space currents and electrode 
currents produced by a sinusoidal potential difference of angular frequency w applied 
between a pair of electrodes at a distance d apart and immersed in a gas discharge. It 
is shown that the ionized gas behaves like an isotropic medium of dielectric constant 
{1 —47ne?/mw?} and zero conductivity only in the limiting case of B=v/wd<1 where 
v is the R.M.S. velocity of the free electrons in the gas. The variation of the equivalent 
dielectric constant and conductivity of the gas with B are sketched. Extensions of Green’s 


reciprocation theorem and of a theorem on induced currents given, by W. Shockley are 
derived in the appendix. 


Thermionic Emission from Oxide Coated Cathodes, by D. A. WRIGHT. 


ABSTRACT. It is shown experimentally that with Ba/Sr oxide cathodes, the saturated 
emission depends on the core material in use, and there is some correlation between total 
coating resistance and saturated pulsed emission. The steady-state D.c. emission is lower 
than the pulsed emission, but the D.c. emission measured immediately on application of the 
anode voltage is similar to the pulsed emission. ‘The difference is due to a decay effect with 
-a time constant in the range 1/1,000-1/10 sec., which is similar in many respects to the decay 
effect studied in an earlier paper on conductivity measurements. It is shown that the appli- 
cation of semiconductor theory can give a consistent account of coating conductivity and 
emission, provided the concentration of barium in stoichiometric excess is about 3 x 1017 
atoms/cm*. ‘This figure has some experimental support. It is shown that the presence of 
an interface layer of the type actually detected should cause decay phenomena with a time 
constant less than 1 psec. and in most cases less than 1/10 usec., so that the decay from pulsed 
to D.C. emission is not due to the capacity effects introduced by the presence of the layer. 

With thoria there is little difference thermionically when studied as a coating on tantalum 
or on tungsten, or as a ceramic tube, and the pulsed emission is about twice the D.c. emission. 
Experiment and theory indicate that thoria is an excess semiconductor containing a stoichio- 
metric excess of about 10'* atoms/cm* of thorium. There is very little sensitivity to 
oxygen, however, at temperatures above 1,850° k 
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